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PREFACE 


The purpose of this contract was to provide data processing and 
analysis in support of two 6SFC programs: tornado detection by analysis 
of radio frequency interference in various frequency bands and sea state 
determination from short pulse radar measurements. During the contract 
period, the contract objectives were completed through various data 
reduction techniques. Strip chart recordings were made of the analog 
tapes collected during the tornado collection exercise. These were used 
for preliminary analysis and as a guide for the selection of the most 
interesting data for digitization and further analysis using digital 
computer data reduction techniques. Computer programs were developed for 
this digitization and analysis. Other computer programs were developed 
and improved for the analysis of the short pulse radar data. Finally, 
computer simulations were developed for the short pulse radar returns 
using Monte Carlo methods on a two-dimensional electromagnetic propagation 
model. Comparisons of the modelled data with measured data indicated a 
high degree of similarity between the two sources. 
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I. Introduction 

The data supplied by GSFC for analysis and data reduction consist 
of 14 track analog tornado tapes and 9 track, 800 cpi digital short 
pulse radar tapes. The tornado tapes were generated by another contractor 
to GSFC under GSFC direction and contain the information shown in 
table I.l.The recordings cover 4 HF and VHF frequency bands with vertical 
and horizontal polarizations and with linear and logarithmic scales. 

In addition the outputs from 2 different lightning stroke detectors 
were recorded. Recording, speeds of 30 ips and 60 ips were used. A data 
base of over 50 tapes was used by Techno-Sciences for the analysis 
and reduction. 

Techno-Sciences provided analog strip charts of all available 
tapes on GSFC supplied equipment in coordination with the contract 

i 

technical officer during the contract period. A complete log of this 
processing is included with this report in table L2.. Because the strip 
chart recorder has only 7 useable channels, only half of the tracks on 
each tape could be done in each pass. However, it was found that 7 
selected tracks had most of the interesting information anyway. 

After the completion of the strip chart recording task, digital 
computer programs were developed for the digitization of selected data 
portions from an analog tape recorder playback. These portions were 
selected under the direction of the GSFC technical officer. The digitized 
tapes then were displayed on various time scales to reveal the detailed 
structure of the lightning waveforms. An analysis program was developed 
to compute interpulse time histograms and amplitude distributions. 


1 


The short pulse radar data was recorded during the JONSWAP exercise 
by 6SFC on 9 track 800 cpi digital tapes using programs developed by 
Techno-Sciences. 23 tapes were generated under a variety of sea conditions 
and equipment configurations. Techno-Sciences developed analysis techniques 
for the recorded data including histogram, spectral density and interpeak 
analysis techniques. 

An important tool in the difficult task of determining sea state 
conditions from short pulse radar data is the use of simulated data. 

Data generated by simulations is highly controlled so that It is 
possible to find and evaluate estimation methods without uncertainty 
as to the true value of the quantity to be estimated, Techno-Sciences 
devloped a 2 dimensional computer simulation during the contract period. 

The modelled data has a high degree of correlation with the measured 
data in appearance. 


TABLE I.l 

TORNADO ANALOG TAPE TRACK ASSIGNMENTS 


Track 

Mode 

Assignment 

1 

Direct 

V-3 Mhz Lin 

2 

Direct 

V-3 Mhz Log 

3 

Direct 

V-30 Mhz Lin 

4 

Direct 

V--30 Mhz Log 

5 

Di rect 

V-VHF 139 Mhz Lin 

6 

Direct 

V-UHF 295, Mhz Lin 

7 

Direct 

H - VHF 139 Mhz Lin 

8 

Di rect 

H - UHF 295 Mhz Lin 

9 

FM 

10 Khz Time Mark or 100 Khz at 60 ips 

10 

FM 

WWVB BCD Time Code 

n 

Di rect 

V-Lightning Stroke Det-Lin 

12 

Di rect 

V-Lightning Stroke Det-Log 

13 

FM 

Taylor Tornado Det-Far 

14 

FM 

Taylor Tornado Det-Near 


TABLE 1.2 





TORNADO 

DETECTION ANALOG TAPE PROCESSING 

1 LOG 


Date Time 

Processed 

Tape 

Designation 

Record 

Speed 

Playback 

Speed 

Strip Chart 
Speed 

Tracks 

Operators 

Comments 

9/17/75 

0920 

8/5/75:1647 

60 ips 

60 ips 

.5 cm/ sec 

1-7 

LD, JB 

Very active 

9/17/75 

1 

i 

1015 

8/5/75:1647 

60 ips 

60 ips 

.5 cm/sec 

8 -14 

LD, JB 

1st minute sporadic 
Tr 14 inactive 
Timing pulses in tr 13 
approx 8 sec period 


TTO 


30 ips 

30 ips 

.5 cm/ sec 



Tr i 3&14 as above, good 
BCD time code on tr 10 

9/\in^ 

1 

8/5//5:1827 

30' i p's 

30 ips 

. 5 cm/sec 

g~rpr 

LD,JB,DL 

Discontinued because of 
doubt about frack^.' 
assignments on tape. 
Printout not kept. 

9/17/75 

1225 

8/5/75:1827 

30 ips 

30 ips j 

.5 cm/sec 

1-7 

LD-DB 

Recorder gain low on 
first half of output. 
Strip recorder out at 
end of tape. 

9/17/75 

1600 

6/29/75: #5 

60 ips 

60 ips 

.5 cm/ sec 

1,3, 5, 6 
8,11,13 

LD,JB 1 

1 

Quiet data. 
Tr 13 blank 

9/17/75 

\ 

6/29/75: #3 

60 ips 

60 ips 

.5 cm/ sec 

1,3, 5, 6 
8,11,13 


Quiet data. 
Tr 13 blank. 

9/17/75 

1645 

8/5/75:1909 

30 ips 

30 ips 

.5 cm/sec i 

1 

1 

1,3, 5, 6 
8,11,13 1 


Very active data. 

9/19/75 

1430 

8/26/75:1630 

#3 

60 ips 

120 ips 

1 

1 cm/ sec 

1 

1,3, 5, 6 
8,11,13 

JB 

1 

Very active, low 
background noise 
Tr 13 blank 

9/19/75 

1455 

8/Z6/75.-1645 

#4 

• 

60 ips 

120 ips 

1 cm/ sec 

1 

1,3, 5, 6 
8,11,13 

1 

u_ 1 

JB 

1 

Very active, low 
background noise 
Tr 13 blank 



TORNADO DETECTION ANALOG PROCESSING LOG 


Date Time Tape Record Playback Strip Chart Tracks ' Operators Comments 
Processed Designation Speed Speed Speed 


9/19/75 

1510 

8/26/75 
# 5 

j 

60 ips 

120 ips 

1 cm/sec 

1,3, 5,6, 
8,11,13 

JB 

Very active Sferics 
Low background noise 
Tr 13 blank 

9/22/75 

1230 

8/26/75 
# 2 

30 ips 

120 ips 

2 cm/ sec 

1,3, 5, 6, 
8,11,13 

JB,LD 

Tr 5 appears to limit 
on all of the 8/26/75 
tapes. Tr 11 has a 
periodic waveform & 
tr 13 has occasional 
spikes, but is other-" 
wise inactive. Trs. 
1,3, 5, 6, 8 all very 
active sferics. 

9/22/75 

1300 

8/26/75 

#1 

60 ips 

120 ips 

1 cm/ sec 

1,3, 5, 6 
8,11,13 

JB,LD 

see above 

9/22/75 

1325 

8/5/75:1853 
# 1 

. 

60 ips 

120 ips 

1 cm/ sec 

1,3, 5, 6, 
8,11,13 i 

JB,LD 

Tr 13 active near end 
Tr 11 has periodic 
pattern. Trs 1,3, 5, 6, 
8 active 

9/22/75 

1342 

8/19/75:1722 
# 1 

60 i ps 

120 ips 

1 cm/ sec 

1,3, 5,6 
8,11,13 

JB 


9/22/75 

9/23/75 

9/23/75 

9/23/75 

1400 

0900 

0935 

1100 

7/8/75:# 1 
7/8/75:# 2 
7/8/75:# 3 
7/8/75:# 4 

' 

. 

60 ips 

120 ips 

1 cm/ sec 

1,3, 5,6 
8,11,13 

JB 

Moderate actvity on 
all channels. Taylor's 
detector (Jr 13) work- 
ing:spikes emitted at 
lightning bursts. 
Lightning stroke det- 
ector (tr 11) puts out 
noise&timing only 



TORNADO DETECTION ANALOG TAPE PROCESSING LOG 

Date Time Tape Record Playback Strip Chart Tracks Operators Comments 
Processed Designation Speed Speed Speed 


9/23/75 

1330 

6/28/75:#! 

30 

ips 

120 

ips 

2 cm/ sec 

1,3, 5,6 

JB 

Not rewound 

9/23/75 

1400 

6/28/75:#2 

60 

ips 

120 

ips 

1 cm/ sec 

8,11,13 

Inactive data 
Tlr 13 not working 

9/24/75 

0925 

6/28/75:#3 

30 

ips 

120 

ips 

2 em/sec 

1,3, 5, 6, 


Inactive data-almost no 


1005 

6/28/75: #4 

60 

ips 

120 

ips 

1 cm/ sec 

8.11,13 

LD 

signal .tape deck noise 
only. Timing marks on 
tr 11, occasional 
spikes on tr 13. 

9/26/75 

0855 

6/25/75:#l 

30 

ips 

120 

ips 

2 cm/ sec 

1,3, 5,6 

JB 

Quiet data. Lightning 


■ 

1656 






8,11,13 


Stroke D§t(tr 11) & 


0930 

:#2 

30 

ips 

120 

ips 

2 cm/ sec 



Taylor's detector(tr 


1030 

:#4 

30 

ips 

120 

ips 

2 cm/ sec 



13) working sporadic- 
ally at best. 


1100 

:#3 

30 

ips 

120 

ips 

. 

2 cm/ sec 



9/26/75 

1325 

6/29/75 :#1 

60 

ips 

120 

ips 

1 cm/ sec 

1,3, 5, 6 


Quiet data. Track 13 









8,11,13 


has a periodic pulse 


1425 

:#2 

30 

ips 

120 

ips 

2 cm/ sec 


JB 

train. 


1500 

:#4 

30 

ips 

120 

ips 

2 cm/ sec 





1520 

:#6 

60 

ips 

120 

ips 

1 cm/ sec 





1545 

:#7 

50 

ips 

120 

ips 

1 cm/sec 




9/26/75 

1620 

7/5/75:#l 

30 

ips 

120 

ips 

2 cm/sec 

1,3, 5,6 


Low to moderate 









3,11,13 

JB 

activity. LS det(tr 11) 


1640 

:#2 

50 

ips 

120 

ips 

1 cm/sec 



works sporadically. Tr 
13 has a periodic pulse 
train 



TORNADO DETECTION ANALOG 


Date Tape Record Playback Strip Chart 

Processed Time Designation Speed Speed Speed 

9/30/75 1000 8/12/75; 60 ips 120 ips 1 cm/sec 

i 1610 #1 

1025 i 8/12/75: 30 ips 120 ips 2 cm/sec 

' ' 1624 m 

9/30/75 1045 8/T8/75: 60 ips 120 ips 1 cm/sec 

1645 #1 

1055 8/18/75: 30 ips 120 ips 2 cm/sec 

1735 #2 i 

i ^ j 

i 9/30/75 1155 9/12/75: 30 ips[l20ips 2 cm/ sec 

‘ 1547 #2 J 

9/30/75 1205 7/6/75 #1 60 ips 120 ips 1 cm/sec 

' i 

9/30/75 1420 9/10/75: 60 ips |120 ips 1 cm/sec 

1645 #1 i 

9/30/75 1435 3/11/75 60 ips |l20 ips 1 cm/ sec 


9/30/75 1500 9/10/75: 30 ips 120 ips 2 cm/sec 

1700 #2 


9/30/74 1510 9/12/75: 60 ips 120 ips 1 cm/sec 

1533 #1 

[530 9/12/75: 60 ips 120 ips 1 cm/sec 

1625 #3 

1550 9/12/75: i 30 ips 120 ips 2 cm/ sec 

J640 #4 

1605 9/12/75: 60 ips 120 ips 1 cm/sec 

1710 #5 i 


PROCESSING LOG 



Tracks Operators Comments 


i 


1 , 3 ,5 )6 J 
8,11,13 

JB. 

JB 

■ 1 

Not gooci data 

1 . 

1 ' 

1 

1 

j 

1 ,3,5,6, 
8,11,13 

JB 

1 

Something on chan. 5 
only. 


JB 

I 

1 ,3,5,6, 
8,11 ,13 

1 

B 

Discopnected from DC 
offset board. Wcrks 
better now. 

1 ,3,5,6, 



Good data 

8,11 ,13 

JB 


1 ,3,5,6, 

1 

U 

B 

1 Good data 

8,11 ,13 



i 


1 ,3,5,6, 
8,11 ,13 

JB 

Very good data 

1 

1 ,3,5,6, 

JB 

i 1 

I Small related pattern j 

8,11,13 


1 nothing on chan. 11 & t 

I 


i 13. i 

1 ,3,5,6, ! 

JB 

* 1 -vw-rrr • - 

• Data on al 1 chan . but i 

8,11,13 i 


! 13. 

1 


i 

j 


! 







































TORNADO 


Date 

Processed 

Time 

Tape 

Designation 

Record 

Speed 

9/30/75 

1615 

7/6/75 #2 

60 ips 

10/7/75 

900 

8/26/75: 
1645 #4 

60 ips 

10/7/75 

920 

9/12/75: 
1547 #2 

30 ips 

10/7/75 

1000 

6/24/75: 
#1 noon 
test 



1140 

6/24/75 
Noon test 
#2 



1100 

6/24/75 

30 ips 



#1 PM 
f 1 i ght 



1120 

6/24/75: 
1729 #2 

30 ips 

10/7/75 

1210 

8/12/75: 
1624 #2 

30 ips 

10/7/75 

1230 

• 

8/18/75: 
1645 #1 
8/18/75: 
1735 #2 

60 ips 
30 ips 

10/21/75 

945 

8/12/75: 
1610 #1 

60 ips 

10/21/75 

1000 

9/12/75: 
1533 #1 

60 ips 






10/21/75 


1020 


7/6/75 #1 


60 ips 


DETECTION ANALOG TAPE PROCESSING LOG 


Playback Strip Chart 
Speed Speed Tracks Operators 


120 ips 1 cm/ sec 


120 ips 1 cm/sec 



120 ips 


120 ips 



120 ips 2 cm/sec 


120 ips 2 cm/ sec 


120 ips 2 cm/ sec 
and 

1 cm/sec 
120 ips 2 cm/ sec 

120 ips 2 cm/ sec 


120 ips 1 cm/ sec 

120 ips 2 cm/ sec 


1 cm/ sec 


1 cm/sec 


1 ,3,5,6, 
8. 11,13 


1 ,3,5,6, 
8,11, 13 


1 ,3,5,6, 
8,11,13.. 



8,11 ,13 



1,3, 5, 6, 
8,11,13 


8,11,13 



Comments 


JB 

Data on all chan, 
but 13. 

JB 

Nothing on chan. 11 

JB 

Nothing on chan. 13 

• 

JB 

Nothing on chan. 13 

JB 

No data 

JB 

Light data 

JB 

Something on 5 & 
11 only. 

JB 

Data on all Chan, 
but 8 and 13. 

JB 

Data on all chan, 
but 8 and 13. Not 
much on 3. 




















































TORNADO DETECTION ANALOG TAPE PROCESSING LOG 




Date Tape 

Processed Time Designation 


Record Playback Strip Chart 

Speed Speed Speed Tracks Operators 


Comments 


10/21/75 1035 7/6/75 #2 


1040 7/6/75 #3 


60 ips 120 ips 1 cm/ sec 1,3, 5, 6, JB 

8,11,13 

30 ips 120 ips 2 cm/sec 


Data on all chan, 
but 8 & 13 

Light data in beginning 
then everything stops 
half way through. 


10/21/75 1300 9/12/75 


varies j as 


,1 cm/sec 1,3, 5, 6, 


5 tapes on one chart 

L. 



II. Analog Strip Chart Recordings 

Table II is a processing log which summarizes the useful analog 
strip chart recording results by Techno-Sciences on GSFC supplied 
equipment. In addition to these results, other recordings were made 
but not kept because of various equipment problems which resulted in 
invalid or useless data displays. Many of the tapes were found to be 
very quiet, containing indiscernable signals or noise of a "white" 
rather than of the desired impulsive nature characterizing lightning 
discharges. Others, however contained considerable amounts of the 
desired impulsive noise across all frequencies. A short segment of 
one of these displays appears in figure II. 1. 

The data was recorded at a real time (one-sided) bandwidth of 
300 khz. As shown in Appendix A, the strip chart recording technique 
used is physically limited to be able to show only about 20% of the 
actual bandwidth. Even this could only be achieved at 1-7/8 ips 
playback speeds and 50 cm/sec strip chart speeds, which would generate 
unacceptably large amounts of output. Hence, the strip chart recordings 
can only be viewed as a "quick look" tool. Selected portions of the 
data must be digitized for processing and display. 

As shown in Appendix B, the amount of digital data generated from 
even a short segment of analog tape is very large. Hence, the segments 
must be carefully selected. The addition of a timing track may be 
necessary in this regard as an aid in registration at some point. 
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figure ll.l. Sample Strip Chart Recording 
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III. Tornado Data Digitization 

Selected portions of tornado data were digitized by Techno-Sciences 
in coordination with GSFC. The Ampex PR-2200 tape recorder used in the 
recording of the data was interfaced to the Interdata Model 5 through 
an A/D converter which was originally developed for data compression. 
Computer programs were written by Techno-Sciences for the digitization 
and recording of the original analog tape data on disk and 9 track 
digital tape. Data recorded on disk is transferred to tape after the 
disk pack is full. 15 bit 2's complement samples were recorded. The disk 
record mode supports a rate of slightly less than 20,000 samples per 
second. The direct tape record mode supports a rate of slightly less 
than 10,000 samples per second. Most data has been recorded directly 
to tape in the interests of conserving operator time and maximizing 
the data throughput on each record cycle. Each disk pack holds approx- 
imately 2.5 megasamples, less than % the storage capacity of a digital 
tape. Thus, in terms of real time, when operating the analog playback 
equipment so that the Nyquist rate of 600,000 samples per second is 
achieved, approximately 20 seconds of real time appear on a digital 
tape in the direct mode whereas only approximately 5 seconds of real 
time can be accomodated on disk record pass. Nonetheless, because of 
the low pass cutoff of the analog tape recorder, it is sometimes des- 
irable to record at the higher rate to disk. 

Appendix C contains a listing of the program developed for the 
digitization of data for disk storage. The program for direct tape 
storage is similar, but is not given. 


IV. Computer Programs for Tornado Digital Data Reduction 

Three major computer programs were developed for the analysis 
and reduction of the digitized tornado data. The goals of the prog- 
ramming effort were to develop the capability of displaying selected 
data portions and analyzing them for their statistical characteristics. 
Thiis required the development of an appropriate printer/plotter routine 
in addition to the formatting, I/O, and analysis routines. Additional 
effort was expended on programming for joint operation on the Interdata 
Models 5 and 8/32. 

The three programs are called SHOALL, REDUCE, and TORNANAL. 

The listings of these programs appear in Appendix D. The plotting 
program, called PLOTIT, appears in Appendix 6‘. Program SHOALL plots 
out the input data, removing non-data portions (noise-only) from the 
plot based on an input threshold and time-window. Other program inputs 
are the number of points to be used in a FFT spectral analysis, the 
number of tape records (one tape record=12800 samples) to be processed, 
the real time sampling rate in kilosamples/second, the time between 
plotted timing lines in milliseconds, and the number of times each 
sample is to be repeated for possible axis expansion. If the repeats 
are set equal to 0, no waveform plotting is done. Similarly, if the 
number of FFT points is set equal to zero, no spectral analysis is 
done. The output plot starts with a header which is entered by the 
operator. This is followed by the waveform plotted as a function of 
time. Everytime a portion of the input is skipped due its detection 
as noise only, the amount of time skipped is printed out. Upon the 
completion of the processing of the desired number of records, or 


upon the detection of tape abnormal status (normally end of file) on 
read, the program prints out the magnitude of the average power spectrum, 
accumulated by the FFT operation over each data portion, and terminates. 

Program REDUCE is for quick-look plotting of the digitized 
tornado data. Inputs are the real time sampling rate and the desired 
output real time plotting rate. The data is filtered by a simple RC 
filter with 3 db point at h the Nyquist rate to reduce aliasing. The 
output plot starts with a header which is input by the operator to 
identify the data, date, etc. and is followed by the plotted data 
with vertical timing lines every 10 msec in real time. The program ter- 
minates upon abnormal tape status (normally end of file). 

Program TORNANAL is used to do statistical analyses on the 
tornado data. An amplitude histogram and a histogram on the time between 
threshold up-crossings is computed. Inputs to the program are the num- 
ber of tape records to be processed, the minimum time interval between 
up-crossings to be counted as a peak event and the up-crossing threshold. 
Outputs are the input parameters and plots of the histograms on a log 
scale. Examining the sample printout in appendix D, it is seen that the 
histogram cdf is approximately linear, indicating an exponential dist- 
ribution for the interpulse arrivals and hence a Poisson peak occurence 
model. The amplitude histogram indicates agreement with a Gaussian amp- 
litude model. 

The reader is referred to Appendix D for listings of the three 
main programs and for sample printouts of each program. 


V . Data Reduction for the Short Pulse Radar Data 

The short pulse radar data is recorded on 800 cpi digital tape. 
Program ECKPDF was developed for the analysis of this data. The listing 
of the program and a sample output appears in Appendix E. The program 
performs amplitude histograms, spectral analyses, and interpulse time 
histograms. In addition, the average pulse shape across the input data 
and the shape of each input pulse can be plotted. These pulses appear 
on a linear scale in time, as’-they are recorded, and on an estimated 
linear scale in range and on an estimated linear scale in antenna 
angle. The necessity for estimation arises from the need for an absolute 

I 

range and antenna angle calibration on the data at some point in the 
recorded time record. This estimated point is found by taking the lead- 
ing edge of the return pulse and assuming that it corresponds to the 
lower 3 db point of the antenna. The range and angle scales are then 
determined from this assumed angle and the known aircraft altitude at 
the time of recording. 

The program inputs are the desired header on the output plot, the 
antenna parameters-pointing angle, beamwidth, and range, the number 
of pulses to be processed, the record format - documentation size, 
data record size, and total record size, and the number of pulse 
points to be used in the reduction. A sample of the program output 
appears with the listing in Appendix E. The input header is followed 
by the record recorded on the tape giving the radar operator information. 
Next appears an amplitude histogram. The non-continuous nature of the 
histogram in the example is due to a misadjusted A/D converter. Next 
appear the averaged pulses in time, range, and angle. Following this 
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is the average pulse power spectrum calculated and accumulated using 
an FFT routine. Finally, the interpulse time histogram is plotted, the 
peak of which can be used as an estimate of the average wave spacing. 
In the example in Appendix E, the last page is a sample of individual 
pulse printouts, on a linear time, range, and angle basis. These are 
optional and would normally appear before the average pulses. 


VI. Short Pulse Radar Simulation 


The problem of sea state estimation from measured short pulse 
radar information is very complex. Particularly complicating the 
problem is the lack of accurate "ground truth" data, i.e. the lack of 
exact knowledge of what an accurate estimate on any particular segment 
of data is. A good way of removing this problem is to develop computer 
simulations where the estimated parameters can be accurately known and 
controlled. 

As a part of the contract effort, Techno-Sciences has developed 
a two-dimensional simulation of the received pulse under a Gaussian 
surface model. A listing of the program appears in Appendix F along 
with a sample printout of several pulses. The surface is allowed to 
consist of any combination of a deterministic sinusoidal component 
and a Gaussian component with a Pierce-Moskowitz spectrum under operator 
control. In each case, the period and amplitude is variable. The random 
portion is generated using an FFT on randomly generated frequency com- 
ponents whose variancsare porportional to the Pierce-Moskowitz spectrum 
at each frequency. The simulation of the pulse return is based on equa- 
tion (4) of Spectrum of Power Scattered by a Short Pulse From a Stoch- 
astic Surface , by David Levine, NASA report X-952-74-299, August 1974. 
The specular reflection points are found by interpolation between sam- 
ples of the derivative of the surface where it goes through zero. 

It has been found that the simulated results agree well with the 


measured data. 


VII. New Technology 

There are no reportable new technology items resulting from the 
work under this contract. The following review activities were per- 
formed to determine any reportable items: 

1. The key technological concepts and ideas studied under the 
contract were identified. These consisted of the methods of analysis 
of tornado and short pulse radar data. The extent to which these ideas 
represented new techniques as versus an application of known tech- 
niques was reviewed. 

2. A review of appropriate published literature to determine 

the uniqueness of the ideas developed under the contract was performed* 

3. A meeting with the technical officer to discuss the results 
of the contract study effort and points (1) and (2) in connection with 
efforts performed at GSFC and under contract with other contractors 
was held. 

As a result of the review activities, it was concluded that there 
were no ideas, discoveries, or improvements or reportable items which 
were first conceived or reduced to practice under the contract. 
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Strip Chart Recorder Display Limitations 

The Varian strip chart recorder used for the processing of the 
analog tape recordings under this contract is limited in usefulness by 
sampling theorem considerations presented in Appendix A. The plotting 
density of the strip chart recorder is 200 points per inch, which by 
the sampling theorem allows for a maximum frequency resolution of 
a waveform of 100 cycles per inch. The maximum strip chart speed is 
50 cm/sec or approximately 20 inches/sec. Thus the maximum frequency 
content of a signal which can be resolved is approximately 2000 hz. 

At the recorded speed, the anlog tapes have frequencies up to 
300 khz • The maximum record speed used was 60 ips. At a minimum 
playback speed of 1-7/8 ips, the highest frequency present is then 
300/32^5; 10 khz or 5 times the strip chart recorder's frequency 
resolving capability. 

Thus digital methods must be used to accurately represent and 
display the analog tape data. From the preceding considerations it 
is seen that a minimum sampling rate of 20,000 samples per second 
must be used. Higher sampling rates are desirable to avoid analog 
tape flutter/wow problems at 1-7/8 ips by running higher playback 
(and hence higher sampling) speeds. 


Appendix B 


Digital Tape Limitations 

Based on the sampling rate considerations presented in 
Appendices A and B, it can be easily shown that great selectivity 
must be exercized in choosing the data to be digitized for display/ 
analysis. Each of the 14 300 khz analog tape tracks requires 
600,000 samples per second in real time or, at a recording speed 
of 60 ips, 10,000 samples per inch of analog tape per channel for a 
total of 140,000 samples per inch for all channels. Assuming 8 bit 
samples and continuous recording (no record gaps), a 2400 ft. 800 bpi 
digital tape can hold only (2400 x 12 x 800)/140,000 = 165 inches 
of analog tape or less than 3 seconds of real time sampled data. Hence 
the channels chosen for sampling and the actual segments of time 
should be carefully chosen for maximum information content. 
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♦ 0PT LA8-FIT0DK 

♦ PR0GRAM T0 INPUT FP0M X'20' A/T< 

♦ 0N SELCH X'Fl‘ AND 0UTPIJT 

♦ THROUGH SVC 1 T0 BRYANT IN 0NE 

♦ CYLINDER CHUNKS ON THE X'FO' SELCH 

♦ SECTOR C0Uii4' AND BUFF BEING FILLED 

♦ ARE DISPLAYED ON D. P. PROGRAM QUITS 

♦ ON ANY INTERRUPT OTHER THAN X'Fl’ DEV 

♦ WITH X'14‘ STATUS (CLOCKED M0DE-EXA) 

4 INCLUDING TTY. NORMAL TERMINATION IS 

♦ THEN V;HEN disk IS FULL. DIRECT 
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EQU 

11 

oooc 
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STOP SELCH 
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READ ONE MORE 
SEND A(BUFF2 > HIGH 


GO SELCH 
SEE IF DONE 

IF SO. GO TO TAPE 

INCR SEC PTR 

SET FOR NEXT BLOCK 

5. DISPLAY 

SET first samps 




• nn'anss!!! 


FI TO DISK UNDER SVC I 11/75 


PAGE 



0172R 





014SR 

§U8r 


3VC 

l.RDD£.< 

READ BLK IN 

014CR 

4310 

0172R 

4230 


LH 

Rl. RDDSK-f'Z 

CK STAT 

0150R 


8N2 

ABORTS 
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TAPMES 

DC 

7. 10. C 

OICOR 

8302 
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LUP 
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NUMSEC 
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PAWS 
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oib5r 
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012AR 

RO 
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OOOB 
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OOOC 
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GOOD 

R14 
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R15 

GOOF 
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0170R 
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0009 

REUDSK 

0132R 

RHIGHl 
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RHIGH2 
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RL0W1 
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RL0W2 
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RST0P 
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SAV 
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Program SHOALL 


PAGE 0001 


fHSSM 

3H&ALL PROG TOPNmDJ PROCESSOR PROGRAM 
ERLST 
IF0RT 

•: #LAB- SHOALL 

C PROGRAM TO READ IN PROM LU 1 THE 
C RECORDED DIGITI2ED TORNADO SAMPLES 
C IN l£.SOO IS BIT SAMPLE BLOCKS AND PLOT 
C THEM ON THE V API AN. 

IMPLICIT INTEGFR*2 (I-N) 

INTEGERfA ISTmI.ISTDEV 
COMPLEX CMFLX 

DIMENSION INF( J 2800 )< IPLOK 3 ). IGD( 21 ), DATT 15 ) 
DIMENSION SPEC< S13 ■. TBU 514 ) 

COMPLEX £( 1024 ■ 

EQUIVALENCE ( IMP( 1 J, DAK 1 ) ) 

DATA C215*I5ER0. I0NE/'327-;e. .0, 1/ 

DATA ICP. IGDFTS, SCALE.-IS-sO, 21. 5. /” 

DATA LU/'j,'' 

GDINC- 1333. .PLOAK I GIFT 3- 1 ) 

UP GRID VALUES 
DO 4 I-l.IGDFTS 
X-I- 1 

IGDf I )»X*GDIf' C + . 5 
'..■RITEC 0- 1 ) 

FORMAK oHL ABEL'S • 

READ (0.2) Di- 7 
FORMAT ( 15A4 ) 

\4'RITEt3. 7) DAT. ICR 
FORMATC IH . ISA^^. A1 ) 

'iVRITE( 0. 21 ) 

FjRMATC 13HNFFT. RECORDS'? ) 

READf 0. 17 ) N, NF ECS 
RECORD COUNTER TO 0 
NRECNT-0 
UP FFT TABLES 

CALL F0URIT< E. M. TBL. IHERO ) 

N2-N/2 
fi21-N2+ 1 

C HERO AVG SPECTRA 
DO 20 I“1.N21 
SPECC I )-0. 

WRITE(0. 14) 

FORMATi; 38HDELAVC MS THPESHC V ). SPACE'’. ). RATE*: KS ) ) 
READ( 0. 15 > DELA.V. THRESH. SPACE. RATE 
FORMAK 15F5. 2 - 

CONVERT TO INTERt'AL PARAMETERS 
RATE IN SAMFS.'SECCSPACE cET^/EEN I'ERTICAL LINES 
IN SECONDS. DELAY Ffi R WINDOW IN SAMPLES. 

THRESHOLD IN INTEGER RANGE 
RATE-RATE* iOCG. 


C SET 


21 

C SET 
C SET 


20 

14 

15 
C 
C 
r 
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SPACE-SFHCE/IOUO. 

IS?AC*SPACE*F A1 E + . 5 
DELAY-DELAt/lQOO. 
idly-delay+rpte*. S 
ITHRSH-( 2. ♦♦ 14 «THRESH/10. 

DT- 1. -RATE 

C GET EXPANSIi^N FACTr.R FOR PLOTTING 
C IRPT-0 MEANS DONT PLOT DATA 
WTTE< 0, \ i » 

10 FQRMAT( SHFPTPEmT ) 

PEhDC 0. IT .[FF TS 

17 FORMAT- ISIS > 

C ADJUST GRID SPAOir-.G FiR REPEATS 
I3 PAC-ISPh<:*IRPT 5 
.'tP.PT-IFFTS 

(.ALL SETGRDC LU, IGD, IGDPT 3, I SPAC ' 

ICNT-0 

SfIP-0. 

ISWIOH-O 

TIME-0. 

• LAST-O. 

NrFT-0 

S CALL SYSI0( Tj. ], I3TAT. ISTDEV. INPT 1 MNPf 12300). 2-0) 

IF( ISTAT. NE. C > GO TO 10 I 

IF( IRPT. NE. 0 . V/RITET 100 TIME 

D0 o III - 1. lIcOO 

TIME-TIMIE-DT 

100 FOPMATC 5HTIME-. F15. 4 ■ 

3RIP-3KIP-DT 

IF( ISWICH. EO. 1 .• GO TO o 

I J-III+IDLY 

IF( IJ. GT. 12SC0 ' G0 TO e 

IF( lABSCINPC IJ • ). LT. ITHRSh ) (,0 T-Zf 6 

ISWICH- 1 

NFFT-0 

SKIPM-SKIP* ICOO. 

TIM-TIPIE* lOOC. 

WRITEC 3. 12 ) SKI PM- TIM 

12 FOPMAK IH -P12.5.21H MSEC SKIPPED- TIME- .F10.2.EH MSEC ; 
GO TO 13 

3 ICNT-ICNT* 1 

IF>: ICNT. LE. lELY ) GO TO 13 
ICNT-0 

LO -I II -IDLY 
Ihl-III-IDLY 

IFTLO. LT. I.OP.IHI.GT. 12300) GO TO 13 
IMA’.<-0 

i'O II 1 1 -LO • 1 HI 
IMX-IABS* IfJF ■ II ' ) 

IllAX-MAKO^ IM; IMO 
IFdMA G G£. ITHFSH) GO TO 13 

11 CONTINUE 
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lSWICH-0 

SK1P«0. 

Gi3 T!3 6 

13 X-INP(III) 

X-X^1C2I5 
NFFT-NFFT* 1 
IF( N. EG. 0 ) <32 T!3 22 
2< NFFT )-Cr'JPL)i( X. 0. ) 

IFCNFFT. LT. N; GS T0 22 
NFFT-0 

CALL FfSUR 1T( r . N. TBL, I!9NE ) 

D0 23 I- 1, N2l 
PE-PEALC Et I ) . 

RI-AIMAG< H( I ; ) 

23 SPECC I )-SPEC( I )+RE+RE+AI*AI 
C LINEARLY INTERP0LATE 

22 IF( IRPT3. EG. C ) 00 T0 6 
D0 18 I- IRFTS 

XX-XLAST + FL0PTC I )*< X-XLAST )/*XRPT 
K-( XX*SCALE*. 5 * 1333. *1.5 
K-liAX0< 0 , K ) 

KK-MIN0( 1333. K ' 

IPL0T( 1 )-KK 

CALL PL0riT< IPL0T. LL. 10NE. I0NE ) 
13 C0NTINUE 

XLA3T-X 

C RESET GRID IF THE 
£ CONTINUE 

NRECNT-NRECNT- 1 
IF< NRECNT. LT. NFECS ) G0 T0 5 
101 IF( N. EG. 0 ) ST0P 

vMo Y_ n 

C FIND AVG Sf=’EC MAX 
D© 24 I - 1, N21 

24 XMAX-AMAXK XMhX. SPECC I ) ) 

C PRINT IT 

D0 25 I-1.N21 
K-SPEC( I 13S3/C(MAX*. 5 
IPL0T< I )-K 
D'3 25 J- 1. 10 

CALL PL0TIT( IPL0T. L. I0NE. ISNE > 


25 

CtiNTINUE 


END 


. U 

EXT 

FUNC 

ISTAT 

INT4 

VAR 

ISTDEV 

INT4 

VAR 

CMPLX 

EXT 

FUNC 

INP 

INT2 

VAR 

IFL0T 

INT2 

VAR 

IGD 

INT2 

VAR 

DAT 

REAL 

VAR 

SPEC 

REAL 

VAR 


OE POOR QUALITY 


> 

J 

1 1 

1 1 1 

l ! 

f 

« 
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T8L 

REAL VAP 




H 

CMPX VAR 




C215 

PEAL VAP 




I2ER® 

I NT 2 VAR 




I0NE 

IHT2 VAR 




ICR 

INT2 VAR 




IGDPTS 

INT2 VAP 




, SCALE 

REAL VAR 




' LU 

INT2 VAR 




GDINC 

REAL VAR 




FL0AT 

EXT FUNC 




FL0AT2 

EXT FUNC 




4 

LABEL 




I 

INT2 VAR 




X 

REAL VAR 




. w 

EXT FUNC 




. r 

EXT FUNC 




i 

LABEL 





EXT FUNC 




2 

LABEL 




•7 

LABEL 




21 

LABEL 




17 

LABEL 




N 

INT2 VAR 




NRECS 

INT2 VAR 




NRECMT 

INT2 VAP 




FOUR IT 

EXT FUNC 




N2 

INT2 VAR 




N21 

INT2 VAR 




20 

LA8FL 




14 

LABEL 




15 

LABEL 




DELAY 

REAL VAR 




THRESH 

REAL VAR 




SPACE 

PEAL VAR 




RATE 

PEAL VAR 




ISPAC 

INT2 VAR 




IDLY 

IN'2 VAR 




I7HRSH 

INT2 VAR 




. R 

EXT FUTX 




DT 

REAL VAR 




\& 

LABEL 




. • IRFT3 

INT2 VAR 




C XRPT 

REAL VAR 




SET '3RD 

EXT FUNC 




ICMT 

INT2 VAR 




SKIP 

PEAL VAR 




IS'UICH 

INT2 VAR 




TIME 

PEAL VAR 




XLAST 

REAL VAR 




NFFT 

INT2 VAR 




* 5 

i 

f 

LABEL 




1 


J) PXGB 





Ob' tooa 
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3YSI0 

EXT FUNC 

lOl 

LABEL 

IRPT 

INT2 VAR 

100 

LABEL 

6 

LABEL 

111 

INT2 VAR 

a 

label 

ij 

INT2 VAR 

IAB3 

EXT FUN3 

IABS2 

EXT FUNC 

SUPII 

PEAL VAR 

TIM 

REAL VAR 

12 

LABEL 

13 

LABEL 

LiZi 

INT2 VAR 

IHI 

INT2 VAR 

I MAX 

INT2 VAR 

1 1 

LABEL 

II 

INT2 VAR 

I MX 

INT2 VAR 

MmXO 

EXT FUNC 

MAX02 

EXT FUNC 

22 

LABEL 

IP 

EXT FUNC 

23 

LABEL 

PE 

real VAR 

PEAL 

EXT FUNC 

AI 

PEAL VAR 

AIMAG 

EXT FUNC 

13 

LABEL 

XX 

REAL VAR 

K 

INT2 VAR 

kk 

INT2 VAR 

MI NO 

EXT FUNC 

MIH02 

EXT FUNC 

PLOTIT 

EXT FUNC 

LL 

INT2 VAR 

. S 

EXT FUNC 

XMAX 

REAL VAR 

24 

LABEL 

AMAXl 

EXT FUNC 

25 

LABEL 

J 

INT2 VAR 

L 

INT2 VAR 

. V 

EXT FUNC 

0000 ERR0RS 
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Program REDUCE - Program Listing 


DEDUCE PRQJ<i Tf0KNrtI>l TIME Ce=MPRE3SIf)N PftStSPftM LOD 
EFL3T 
iFiTiPT 

C TLfTB- REDUCE 

C PRaaRnti Tld READ IfJ TCRNAC® DATA 
■; RC RILTER WITH Tf-E 3 D8 FaiNT AT 
C M'2 THE NtOUIST RATE F'JiR THE OUTPUT 
C PLOTTING SAMPLING PATE 

THE AiJTPlJT RATE IS AN INPUT ^APInBLE 
: :m ksamps/'SEc 
■; PLOTTING CONTINUES T0 EOF 
•: LDD 12'S.''^S MOD FUP a,C52 Z/S.'76 
IMiPLICir INTEGLR»2 'I-N) 

INTEGER*4 ISTAT.ISTDEV 

DIMENSION INF< 12S00 j. IPLOTf 10 ). 1 IDf 21 ) 
DII'IENSION DAT< 15 ). LASTC IQ ; . 

DATA IGDPTS.-^L' 

DATA LUL lONE/T?, I. ' 

GDINC- 1353. /FLOAT'; IGDPTS- i 
C SET UP GRID VAL.UE5 
DU 4 I • 1 ^ I GDF T S 


4 

1 

2 


f 


2 1 
15 



SO 


IJ 1 


<-I- 1 

I'SDT I )- <*'5DITC*. S 
LTrlTE^ 0* 1 ) 

F ORMAT! oHLABELT ) 
r<EAD( 0. 2 ) DAT 
FORMAT! 15A4 ) 

WRITE! 3.. 7) DPT 
FORMAT! IH .15A-) 

‘.4PITE< 0. 21 

FORMAT! 25HFATE IN, RATE OUT.RhNGE'? 
PEAE" 0, 15 ' RPTIN, RAT OUT, RANGE 
FORMAT! 15F5. I ■ 

'.riTEf 3,77 .' FATIN.RATOUT.FAN'iE 
FORr-’AT'' SHF ATE IN-. F3. 3, 4H i<S, • '■•|if?AT 
1 lOHP-P SCrtLE-. F5. 3. 3H V. '■ 

SCALE- 20. ♦ 133 3. - •RANGE 
ALPDEL-3. 14 l = 327*RAT'0UT/-PATIN 
RI-, 0-EXP' -ALFIEL ) 

I3END-RATIN,'RATn)lJT *. S 
ISP AC - 10. ♦RATOUT 5 
SEC 'GRID LINES 
W'PITEf 3, SO > 

FORMAT! 21H. 01 TECOND 'GRID I.IITE3 • 
CALL SETGRD! LU. IGD- IGDPTS, ISPAC ' 
rCNT -0 
IREC-O 


) 

E OUT-. 


2 .' 7 ® 


F5. 3, ah 


T'NURM-' 2. *♦ IE '• SORT'. 1. -FHi5*RH0 • 

11 CALL SySI0! 72 . ; . I’VTAT . I TTDEv' . I'NF' I 1, IW^( 1^:300 • • 2.. 
IF! ISTAT, EQ. C ) 'jO to 3 


K 
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c 1 


^.MTtCO.S) 

FP'Rr-WT( 3H5 ThTUS -.14) 

STitP 

1REC-IREC»1 
URIT£(3,8.'' If- EC 
F‘bRr -.(^TC 7HREC2RL .15) 

Dy 6 III- 1. IcSOO 
X-INP( III ) 

V-RHS-Y+X.-'XNeRit 
ICNT-ICNT+ 1 

IF*; ICNT. LT. I5END ) Qia 4 
ICNT-0 

K-y*SCPLE*^95. 5 
P0INT.lREPEPT.CRID 0N 

CPLL PL0TITC ► . L. I'5N£. I5INE ) 
ISET-ISET + 'l 
CONTINUE 
Q0 T0 11 
END 


. u 

EXT FUNC 

ISTAT 

I NT 4 VAR 

ISTDEV 

INT4 VAR 

INP 

INT2 VAR 

IPL'ST 

INT2 VAR 

IGD 

I NT 2 VAR 

DhT 

PEAL VAR 

LAST 

INT2 VAR 

I*?DPTS 

I NT 2 VAR 

LU 

I NT 2 VAR 

ICME 

INT2 VAR 

GDINC 

REAL VAR 

FLSAT 

EXT FUNC 

FL£*AT2 

EXT FUNC 

4 

LABEL 

I 

INT2 VAR 

X 

REAL VAR 

. U 

EXT FUNC 

. Y 

EXT FUNC 

1 

LABEL 

'I'H 

EXT FUNC 

2 

LABEL 

7 

LABEL 

21 

LABEL 

IS 

LABEL 

RAT IN 

REAL VAR 

RATIMJT 

FEAL VAR 

RANGE 

REAL VAR 

77 

LABEL 

SCALE 

PEAL VAR 

ALPDEL 

REAL VAR 

RH'3 

REAL VAR 

EXP 

EXT FUNC 


( 


V 
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00i?2 
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I SEND 

INT^ ¥AP 

ISP AC 

INT2 VhP 

SJD 

LABEL 

SETdkD 

EXT FUNC 

1C^4T 

I NTS VAR 

IREC 

I NT 2 VAP 

t 

REAL VAR 

1 

PEAL VAR 

. p 

EXT FUNC 

SU:F.r 

EXT FUNC 

1 1 

label 

SY3I3 

EXT FUNC 


LABEL 

s 

LABEL 

.*s 

EXT FUNC 

8 

LABEL 

s 

LABEL 

III 

ir«T2 VAR 

K 

INT2 VAP 

PLSTIT 

EXT FUNC 

L 

INT2 VAR 

ISET 

INT2 VAP 

. ¥ 

EXT FUNC 

0000 EPP'iiRS 



;iur 










CtCiCj t\i or< r<iv •- o o ooooor> 


IrtSSH 

T0RUANAL FR0G T0RNAD'i 
EPLST 
tFORT 

PRfiJGRAfI T0 READ IT 
PECQiPDED DIGITI5ED 
IN l£300 15 BIT SAi 
AND ANALYSE THEM F 
TIME INTERVAL DIET,- 

Iltf’LICIT INTEi.: 
lNTEGEP+4 IPPRI 
DIHENSIijN INP< 
DIMENSION IQPD' 
EQUIVALENCE CIi 
DATA LU/3/ 

DATA IRPTS/5/ 
DATA IGDPTS/c/ 
DATA I0NE, IMl 
DATA ITMMAX, lAi 
AMPLITUDE Normal I S.i 


Pjrogrdm TORNANAL - Program Listing 

DATA ANAL.Y3IS FRSG TSI 4/76 


N0RM- 32763/1 AMtf MX 
SET UP GRID POINTS 
D0 15 I- 1. IGIP 
> IGRD( I )-FLS)AT( 

CALL SETGRDTLU 
V/PITE( 0. 1 ) 

FSRMAH 27HREC0 
F:EAD(0.20) NPE< 

3 F'SRtlATC 215. FE . 

CraNVERT TO INTEGE:: rH:7E.'HC'LD 
ITHR3H»( 2. ** 14 I -THRESH. 10. 

0R - 10. V PANCE 
CLEi-iP PI1F ARRAYS 

L '4 2 I-l.ITMMA 
T iriPMFi; I )-G. 

IJ 3 I-l.IAMFM 
Al ipf ir'! I “G. 

PEAK ORIGIN T0._ 

PEAK E'/ENT CM it ALSO 
IPK0R';»O 
I EVENT*-: 

LAST ONE IN T0 
LASTIN«0 

DO 4 II ;t * I- : r tcs 

CALL SV.-IVjt l-r.TMl 


SET 

AND 


iET 


FROM LU 1 THE 
TI3RNAD0 SAMPLES 
:PLE BL0CKS 
R ATFLITUDE «. 

IEUTI0NS 

P.#2 CI-N) 

IK 6 ) 

2300 ). TIMPflF*. 1024 ). AMPI-F< 1024 ) 
10 ). IDAT'. 5 4 IFARMC 12 ). LAST( 5 3 
ARM< 1 3. IPhRMK 1 3 3 


nW 0 /l, - 1 . 2 / 
IFMX/1024. 1024/ 
R F0R PMF 


- 1 3# 1339. .FL‘3ATC IGDPTS- I 3 
IGRD. IGDPTS. 0 3 

DS. DUR. THRESH?! 215. F5 3 3 
3. IDUR. THRESH 


C EXIT 


IPARMi. 2 I E 

DO 6 III* 1. IcCr 
I PKARG* I FkOF'i ■l■ 
I£v■ENT.IE■''ENT* 
K*( INP( III )- 16’ 


0 ) G a 


1. :Ni-> 

Tj 3 


1 ). 25600. 0 3 


34 VNGPM 1 


I 


ORIGINAX; PAGE 
OP POOR 


QUALm 


I 


X 


hmpmfc k )» ftMPrr 
IF( INP( III ). LT.J 
IF(.LrtSTIN. 'jr. I 
C IF N0T, THIS HN tP 
IF( IPK0RG. LT, I. 
C SKIP SH0RT INTERVhI 
lEVENT-IEVENT 
lEVENT-MINGC IE' 
C N0KMHL12ED Tl3 EVENi 
TIMPMFC I EVENT ■ 
lEVENT-Q 
1* IPK0R(5-O 

6 LASTIN»INP< III 

4 CONTINUE 
'50 T0 12 

C A&N0RMHU END 

5 V.'RITEC 0. 13 ) II 

13 FiaRt-WU 1OH0N RJ 

12 CaNTINUE 

C NWRMALI5ES.PL:3T PMK 
SUM-0. 

m S I- 1* ITMHh 
3 SUM-SUM+TIMPMF 

URITEC 3. 13) NR3 
IS F'3RMHT( IH .1* ^ 

1 F12. 0. 21H TinJ 

i;rite( 3. is ) 

19 F jRMAK 20HL0G 1 

ViRITEC 3. 21 ) 

21 FaRMrtK 2H-5. >3 

CDF- 1. 

LH3U 2 )-0 
LH3T( 1 )-0 
D0 9 I-1>ITMT!m 
TIr!FMF( I )-Ti:r; 
II'hK 1 )-230. f' 
'.DF-CDF-TIMPriF 
IDHTt 2 )-230. r( 
013 S J-l.IRPTS 
CmLL PL'STITklDi 
3 C'jNTINUE 

SUM-0. 

D'3 10 I-l-lMfiPI, 
10 SUM- SUIT* PiiiF MFC 

ClF- 1. 

LtiSTC 1 ■‘-0 
LHSTC 2 )-0 
DO 11 I-l,IHt!P' 
HMPMFC I )-riMFtiF 
IDhTC 1 )«2£0. tt. t| 
CDF-CDF-HMFMF'; 

I DHT < 2 )-230. ♦ C 


Pflr.F 0002 


K )+ 1. 

ITHR5H) GO TO <5 

riRSH) ca TO 6 

CROSSING 

UR ) Ga TO 13 

SCPHRT OF SAME EVENT) 

DUR 

ENT. ITMMHX ) 

LENGTHS 

TIFFMFC lEVENT >r 1. 


I. IPHRMC 2 ) 

CORD .I5> 9H.3THTUS- .1-5) 


I ) 

';S, THRESH, SUM 

OH RE'.ORDS ,.11H THRESHOLD-. FS. 3, 
INTERVALS FOUND ) 

OF PR'O&ABILnv ) 

. 2H-4. 33X. 2H-3. 33X. 2H-2. 33;<. 2H- 1 ) 


Tc I )/sun 

LOG I'jv AM:Y.< IC . 0000 1, TIMPMF' I ) ) ''■* = . ' + . S 
UOG ICC AMHX U . 0000 1. CDF ) )*5. !'+. S 


44 


T . LAST .IT WfD. I i’-T. 
IX 


TX 

I )/3JM 

LOGIO'. orw-f 1C . 00001. Af'TMFC I 0 + 5. X*. 5 
LUGIOC AMH..1C . 00001. CDF ) )+5. >+. 5 


ORlGtNADTAG®? 
OF fOORQ^^' 


Dia 11 J»1.IRFT 
CALL PL'jTITdl*! 


1 1 

CONTINUE 

ST0P 

£|MD 

. U 

EXT KUNC 

IPARlil 

INT4 VAR 

IMP 

INT£ VAR 

TIMPMF 

REAL VAR 

AMP Mr 

REAL VAR 

X 'sPD 

INT2 VAR 

IDAT 

INT£ VAR 

IPARM 

I NT 2 VAR 

LAST 

INT2 VAR 

LU 

I NT 2 VAR 

IRFTS 

INT2 VAR 

I-5LPTS 

INT2 VAR 

ICiUE 

INT2 VAR 

im 

INT2 VAR 

iTv;a 

INT2 VAR 

IThlttiX 

INT2 VAR 

lAMPhX 

INT2 VAR 

Mi-iRM 

INT2 VAR 

IS 

LABEL 

I 

I NT 2 VAR 

PLEAT 

EXT FUNC 

FL0hT£ 

EXT FUNC 

. T 

EXT FUNC 

SET'aRD 

EXT FUNC 

1 

LABEL 

C'H 

EXT FUNC 

20 

LABEL 

NPCCS 

INT2 VAR 

IDUR 

INT2 VAR 

THRESH 

REAL VAR 

ITHRSH 

INT2 VAR 

. R 

EXT FUNC 

£ 

LABEL 

w 

LABEL 

IPKERG 

INT2 VAR 

lEVENT 

I NT 2 VAR 

LAST IN 

INT2 VAR 

4 

label 

IIII 

I NT 2 VAR 

SYSI3 

EXT FUNC 

S 

LABEL 

o 

LABEL 

III 

I NT 2 VAR 

j,; 

INT2 VAR 

I te 

LABEL 

MI NO 

EXT FUNC 

MI Nut 

EXT FUNC 


ooo: 


45 


r. LAST, ITV/<5, IfjNE } 
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la 

LABEL 

13 

LABEL 

SUM 

REAL VAR 

3 

LABEL 

13 

LABEL 

13 

label 

ai 

LABEL 

CE-F 

PEAL VAR 

3 

LABEL 

ALSGIO 

EXT FUNC 

AMAX 1 

EXT FUNC 

J 

INT2 VAR 

PL0TIT 

EXT FUNC 

10 

LABEL 

1 1 

LABEL 

. 3 

EXT FUNC 

. V 

EXT FUNC 

0000 errors: fo 


LEVEL 1 P03-Q0 




I 


Probability Mass Function 
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Appendix E 
Program ECKPDF 

Program Listing and Sample Printout 


Program ECKPDF - Program Listing 

f ASSM 

ECKPDF PROG PDF DATA ANALYSIS PROGRAM E/7£ 

ERLST 

fFORT 

IMPLICIT IHTEGER)K2 (I-M> 

INTEGER5iC4 IPARM(G) 

INTEGER)i<4 IHEDCK , IHED < 40 ) 

DIMENSION INNC 1030) , IHP(204S) , DATC4C) , ISTaT(£) 
DIMENSION IPLOT C 1400) , LAST( 10) 

DIMENSION PULSE(20^S) , VAR(2043) 

DIMENSION PDF(25S) , IGRDC10) 

DIMENSION INEU<1400) 

DIMENSION TBLCB12) 

COMPLEX 2(512),CMPLX 
DIMENSION PAVC512) 

EQUIVALENCE ( Z ( 1 ), PULSE ( 1 )),< PAV ( 1 ), VAR ( 1 ) ) 

EQUIVALENCE C INN C 1 ) , IHEDCK) , ( IPARM( 1 ) , ISTAT( 1 ) ) 

DATA lONE, IZERO , IS/ITUCL 1 3 , 1 25G/ 1 , 0 , 8 , ‘E*, 3 , 25£ / 

DATA ISPAC, IRPT , MAX'PLS/1 2S , S , 2SG/ 

DATA IBIG., IAMP/30000, 2/ 

DATA DELTATZ.01/ 

C DELTAT=SAMPLING PERIOD ( MSEC ) 

DATA NFFT/-258/ 

DATA LU, YE3/3, IHY/ 

DATA DTOR, DLT300/ . 01745329, 3 . / 

) DTCR=DEGREES TO RADI AHS . DLT300=300XSAMPLING PERICDC MICSEC ) 
DATA IS, IMl/6 , -1/ 

DATA IGRDCi) , IGRDC2) , IGRD(3) , IGRD(4), IGRDCB) , IGRDCS)/ 

1 0,279,559,83.9,1119,1399/ 

C GET IHEDER IN 

CALL SYSIOdPARM, 72, 1, 1HED< 1 ) , 80, O) 

CALL F0UR1TC2,NFFT,TBL, IZERO) 

C SET UP FFT TABLES 

NFFT21=NFFT/2+l 

URITEC0,40) 

40 FORMAT <12HY FOR PULSES) 

READ (0,50) PRPULS 
50 FORMAT CAD 

WRITE (0,4) 

4 FORMAT (SH LABEL) 

READ(0,5) DAT 

5 FORMAT (40A4) 

WRITE (0, 30) 

30 FORMAT (2SH ANGLE , BMWI DTH , ALTUD i M ), FE . D 

READ (0, 3D ANGLE, BEAM, ALTUD 
C CONVERT TO RADIANS FOR LATER TRIG 
ANGRAD=ANGLE^DTOR 
BMRAD = BEAM.-«DTOR 
C ONE SIDED BEAMUIDTH 
BM2=BMR AD/2, 

C ADJUST MIN ANGLE FOR NADIR OR NOT 
ANGMIN = AMA,X1 (0 . , ANGRAD-BM2) 
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51 

C ALSO MAX 

ANGMAX>ANGRAD+BM2 

C DXR= X INCREMENT ALONG THE GROUND LJHEN ADJUSTING PULSE BASE 
C CHOSEN SO THAT NO. OF GROUND SAMPS==NO. OF TIME SAMPS 
C FROM NADIR TO ANGMAX 

DXR = SIN( ANGMA.X ):l;DLT300.'-<a . *ALTUDX( 1 . -COS C ANGMAX ) > ) 

C NOU DO SAME FOR UNIFORM ANGLES 

DANG= DLT300:l: ANGMAXr < 2 . T. ALTUD^t^ < : . /COS ( ANGMAX ) - 1 . ) ;> 

31 FORMAT C SFS . 1 ;• 

C CLEAR PDF ARRAY 

DO 43 1=1,356 
42 PDFtI)=0. 

C CLEAR AVG PULSE ARRAY 
DO 60 1=1,3048 
VARC I )=0, 

60 PUL5E(n=0. 

URITEC0, 1) 

1 FORMAT C7H FILE=?) 

READ C0,SY IFILE 
3 F0RMATC8I5) 

URITEC0, 7) 

7 F0RMATC34HBL0CKS,BLKSIZE,D0CSI7E,LRECL,L0,HI ) 

READC0,3:t IBLKS, IREOS, IDOCS, IDaTS, LC, IHI 

IDOCSl =l-f IDOCS/3 

IREC3=IRECS/3 

IDATS£=IDATS/3 

LO = MAX0>: 1, LO) 

IHI=MIN0( IHI , IDATS) 

IHI=MAX0CNFFT+LO-1 , IHI ) 

IHIGH=IHI-LO+l 
IHIGH=MIN0C IHIGH, 45G) 

IHIGH3=3XIHIGH 
0 ALLOU FOR 3 TRACES 
IHIGHU = IHIGH:i<3 
SAMPS= IHIGH 
DUR=DELTAT.TSAMPS 
I BEG = IDOCSl 
KTOT=G 

DO 3 IIII=i, IBLKS 

IF( IBEG, EQ .IDOCSl )CALL SY3I0 C IP ARM , 73 , 1, INN( 1) ,2060,0) 

IFC ISTAT(2) , EQ .0) GO TO 128 
URITE(0,S0) ISTATCa) 

S0 FORMAT C7HSTATUS= ,16) • 

GO TO as 

123 IFC IHEDGK. , NE , IHEDCl J ) GO T0S2S ' 

URITEC0, Si ) 

81 FORMAT CISHHEADER FOUND) 

GO TO 38 

C TERMINATE ON LAST IHEDER FOUND TOO 

338 CALL BYTOHUC INNC IBEG) , IDATS, INPC 1 ). 18 ) 

IBEG=IBEG+IDATS3 

IF C IBEG , GT , IREC3) I,BEG = IDOCSl 





M 
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F0RMATC2EI4) 

KT0T=KT0T+1 
DO 17 K=LO, IHI 
I=INPCK)+1 
X=I-1 

PiJLSECK) =PULSE(K)+X 

varck)=xxx+'v‘arck:i 
PDF( I)=PDF( I )+l . 

FORMAT (4E20 . 9) 

CONTINUE 
SUM=0 . 

DO S 1=1,256 
SUM=SUM+FDF( I) 

XMA>< = 0. 

DO 10 1 = 1/256 

PDFt;i^=PDFCI)/SUM 

XM AK = AM AX 1 C XM AX , PDF Cl ) ) 

URITECS/S) DAT 
FORMAT (IH ,40A4) 

' OUT TAPE HEADER STUFF 
URITE(3,6) IHED 

UR I TEC 3, 33 I F ILE . ANGLE , DUR , SAMPS , ALT UD 

FORMATCSH FILE , I 3 , 5X , F6 . 1 , SH DEGREES , 5X , F6 . c , 7H MICSEC, 
15X,F7,2,SH PTS. IN /iOH , ALTITUDE= , FS.0) 

URITEC3,20) KTOT,LO,,IHl 

FORMAT C IK .. I6,7H BLOCKS , 8H , LI M I TS = , 2 16 ) 

URITEC3,213 

FORMAT C20HLOG OF PROBAEILITV 5 
URITEC3,22) 

FORMAT C2K-5, 33X , 2H-4 , 33X , 2H-3 , 33X , 2P-2 , 33X , 2P- 1 , 33X , 1H0) 
CALL SETGRDCLU, IGRD, 16, ISPAC) 

DO 12 1=1,256 

IDAT=2S0 . XC ALOG10CAMAX1 ( . 00001 , PDF/ I ) ) 3+5 , )+ . 5 
DO 12 J=i , IRPT 

CALL PLOTITCIDAT, IDAT, lONE, IMl j 
UR I TEC 3, 63 3 

FORMAT C20H1 AVERAGE PULSE, SIGMA ) 

' AVG PULSE,, PEAK 
PK = 0, 

XTOT=KTOT 

DO 61 I=LO, IHI 

PIJLSEC I ) =PLILSE( I )/XTOT 

VARCI 3 =VARC I )'/XTOT-PULSEC I3XX2 

IFCPK.GT.PULSECD) GO TO 61 

PK=PULSECI) 

KSTART=I-4 

CONTINUE 

D£LTAR=ALTUDX( 1 . /COS C ANGRAD+BM2 ) - 1 . /COS C ANGMIN ) 3 

DELTAS =2 , XDELTARXDLT300 

KEND = DELTAS-I- , 5 

KEND=XENB+KSTART 

yRITEC3/151) KSTART.KEND 
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151 F0RMATC7H5TART =,I6,7H, END =,];S) 

C FIRST FIND ZERO 
PD0=0, 

DO 51 1=1,356 
PD0=PD0+FDF( I ) 

IFCPD0 . LT , .0801 D GO TO 51 
KZERO=I-i 
GO TO 53 
^ 51 CONTINUE 

: S3 CONTINUE 

: C NOU FIND ADJUSTED AVERAGE PULSE 
XZERO=KZERO 
PK=PK-XZERO 
DO 64 K=LO, IHI 
PULSECK ) =PULSECK.)-XZERO 
64 CONTINUE 

C SET GRID FOR AVG PULSE 

CALL SETGRD(LU, IGRD . 16, IBIG) 

SCALE=13S9./PK 

DO 63 K=LO,IHI 

IPLOTC 1 ) =SCALE'4CPULSECK)+-. 5 

IPLOT(2)=SCALEXSQRT(AMAK1(0, , VAR(K)) )-l-.5 

IFCK.NE,LO:) GO TO 72 

LASTC l)=IPLOTCl) 

LAST C E ) = IPLOT C 3 ) 

- 73 CONTINUE 

: 0 FORCE GRID AT KSTART 

IF(K.EQ.KSTART.OR.K.EQ.KEND) CALL SETGRD ( LU , I GRD , 16 . IBI G 1 
DO S3 J=1 , IRPT 

CALL PLOTITC IPLOT, LAST, ITUO, TONE) 

S3 CONTINUE 

C ADJUST TINE BASE 
C FOR UNIFORM SAMPLING ON GROUND 
C GET START SAMPLE VALUE ON INPUT (MAV EE -) 

DRSTRT=ALTUDXC 1 . /C03C ANGMIN)-! , ) 

K0=KSTART-IFIX(3 . XDRSTRT/DLT300+ . 5) 

C NORMALIZED RANGE VALUES USED 
C DTOT=TOTAL SURFACE COVERED 

DTOT = DXR;t:FLOAT( IHIGH.-)ii?ALTUD 
URITEC3, 35)ALTUD,DT0T 
LAST(1)=0 
KSET=0 

DO 4S K=l, IHIGH 
XR=DXRXFLOAT(K-l) 
i , RNS=1 . +XRXXR 

i . , DR=ALTUDXCSQRT(RNS) -1 . ) 

BS=2 . XDR/DLT300 
KD = DS 

: C XR=NORMALIZED X , RNS = NORMALI ZED RANGE , DR= RANGE CHANGE 

: C KD = TRIJNCATED SAMPLE CHANGE 

KADJ=K0+KD 

IFCKADJ-. NE . KSTART . AND , KADJ . NE . KENDIGO TO 5S 


nL 







r 


IFCKSET . EQ , 1 . AND . KADJ . EQ . KSTARV ) QO TO 49 
KSET=1 


CALL SETGRDCLU, IGRD, 16, IBIG) 

SS IFCKADJ.GT.LO.AND.KADJ.LT, IHI) GO TC 49 


IPLOT(1)=0 
GO TO 57 

49 I PLOT ( 1 ) = SOALE'^ ( PULSE < K AD J ) + ( PULSE ( K AD J+ 1) - 

1 PULSE (KADJ) ):S<DS-FLOAT(KD) J)+ 5 
57 DO 4S 1=1 , IRPT 

CALL PLOTIT<IPLOT,LAST, lONE, IONE3 
4S CONTINUE 

C HOU PUT OUT UNIFORM ANGLE PULSES DANG=DELTA ANGLE 
KSET=0 
LAST ( 1 ) =0 

ANGTOT = DANG'4CFLO AT ( THIGH) 

ANGTOT = iS0.;{CANGTOT/3. 1415927 
URITE(3,GG) ANGTOT 

SS FORMAT CSH 10 TO ,F9.2,SH DEGREES?) 

DO G7 K=1,IHIGH 
AMG=DANG:<FL0AT(K-1 ) 

DR = ALTUD-li C 1 . /COS i ANG ) -1 . ) 

DS = 2 . ;JfDR/DLT300 
KD = DS 

C ANG=AMGLE FROM NADIR , DR=RAHGE CHANGE FRCM NADIR, 

C KD= TRUNCATED SAMPLE CHANGE 
KADJ = K0 + K.D 

IF(KADJ.NE,KSTART, AND.KADJ.NE.KEND) GO TO 68 
IFCKSET.EQ.l . AND.KADJ.EQ .KSTART) GO TO 69 
KSET=1 

CALL SETGRDCLU. IGRD, 16, IBIG) 

C PUT VERTICAL LINE HERE 

68 IF(KADJ . GT , LO . AND , KADJ . LT . IHI )GO TO 69 
IPLOT(l)=e 

GO TO 79 

69 IPLOT(l)=SCALE;<<(PULSE(KADJ) + (PtJLSE(KADJ + l)- 
1 PULSECKADJ) ))«(CDS-FLOAT(KD) ) )+ . 5 


79 DO 67 1=1, IRPT 

CALL FLOTITC IPLOT, LAST, lONE, lONE) 
67 CONTINUE 

C SET UP GRID ^ 

CALL SETGRDCLU, IGRD, IZERO , MAXPLS ) 
URITEC0, 65) 


ORMAT C25HPOSITION TAPE FOR PULSES) 


PAUSE 


C PUT LINE AT PULSE START, END 
IPLOTC1)=0 
, IPL0T(2)=IHTGH-1 
' IPLOT C3) =IHIGH2-1 
C GET MAX RANGE ADJUST FACTOR 
C R'-lc:i?4 FACTOR 

URITE(3,E5) ALTUD,DTOT 

25 FORMATClcHlALTITUDE = ,F7.0,2H M , 221- DURATION FROM NADIR 
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C CLEAR AVG SPEC ARRAV 
DO 89 1=1, 5ia 
39 PAV(I)=0. 

C CLEAR PROB ARRAY 
DO 77 1=1,356 
77 PDF(I)=0. 

MAXAMP = aS6:*;iANP 
MAXAMl =MAXAI1P-2 
TIMINC = ITf1INC 
IBEG=IDOCSl 
DO 70 IIII=1 , IBLKS 

IFCIBEG.EQ. IDOCSDCALL SVSIO < IPARM , 72 , 1, INN ( 1 ) , 2OS0 , 0 ) 

IFCISTAT(2) , NE . 0)GO TO 78 

CALL BYTOHUCINNC IBEG) , IDATS, INP( 1 ) , IS) 

IBEG=IBEG+IDATS2 

IF ( IBEG . GT. IREC2) IBEG=ID0C51 

MAX = 8 

MIN=MAXAMP 

DO 71 K=LO, IHI 

I = MAX0t;0, INPCiO-KZERO):t:lAMP 

J=K-LD+1 

INELJCJ) = I 

MIN=MIN0(MIN, I) 

HAX=MAX0CMAX, I ) 

71 INPCK)=I 

C MAX=MAX VALUE OF ADJUSTED DATA , MIN^MINIMUM 
DO 24 K=1,IHIGH 
XR=DXRXFLOAT(K-l ) 

C NEU X VALUE 

ALF=ATAN(XR) 

X=1 , 39X(ALF-ANGRAD)/BM2 
SX=ABS(SIN(X) ) 

IF(SX,GT.X/2. ) GO TO 26- 
ANTFAC=i, 

GO TO 27 

26 ANTFAC= (X/SX)XX4 

C GET RANGE NORMALIZED SQUARED 

27 RNS=1 , +XRXXR 

C GET RANGE AND ANTENNA ADJUSTMENT FACTOR 
RADJ=RNSXRNS 

C DR= ACTUAL RANGE TO GROUND POINT - ALTITUDE 

dr = altudxi:sg!RT(rns)-i . ) 

DS=H . fBR/DLT300 

C TRUNCATEDt" SAMPLE TIMECMAV BE -) 

KD = DS 

KADJ=K0+KD 

C INDEX OVER FOR TIME PULSE 
J=K+IHIGH 

IF ( KAD J . GE , LO . AMD . KAD J . LT . IHI ) GO TO 91 
C IF -, FILL IN UITH 0 
INEUC J ) =0 
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GO TO 24 

DX=INP(KADJ+1)-INP(KADJ) 

7 GV INTERPOLATE 

f ■ INELI ( J J =RAD JX ( FLOAT C INP ( KAD J ) ) -DX:*: ( DS--FLOAT ( KD ) ) j + . 5 

MIN=MIN0(MIN, INEUC J) ) 

MAX = M AX0 C MAX , I NEU C J 1 } 

24 CONTINUE 
C SET FIRST PEAK INDICATOR 
NOPEAK =0 
DO 3i7 1 = 1 , NFFT 
J=I+IHIGH 

IF(4*INEUC.J) . LT . MAX) GO TO 95 
C ABOVE THRESH, POSSIBLE PEAK 

IFdNEUC J) ,GT . INEU( J+1 ) ) GO TO 95 
C DOWNHILL, DOESNT COUNT 

IF( INEldC J+1 ) , LE , INEUC J+E) ) GO TO 95 
C STILL UP HILL 
IPKN=J+1 

IF CNOPEAK . EQ . 0) GQ TO 94 
C FIRST PEAK 

K = IPK.N-NOPEAK 
K=MAX0( lONE, K) 

K = MIN0 C IE'5S, K ) 

PDFCK) =PDFCK)+1 . 

94 NOPEAK=IFKN 

95 X=INEU(J) 

X=SQRT( AMAXl (0 . , X) ) 

97 ZCI ) =CMPLX(X, 0 . ) 

CALL FOURITCZ, NFFT, TBL, lONE) 

DO 3S r=l,NFFTSl 

98 PAVC I > =CAESCZCI ) )XXa+PAV( I) 

IF(PRPULS . NE, VES) GO TO 70 
DO 84 K=1,IHIGH 
ANG=DANGXFL0AT(K-1 ) 

DR=ALTUDX(1 , /COS( ANG)-1 . ) 

DS=E . XDR/DLT300 

KD = DS 

KADJ=K0+KD 

J=K+IHIGH2 

C PUT THIS LINEAR ANGLE PLOT 3D IN LINE 

IFCKADJ . GE , LO , AND . KADJ . LT , IHDGO TO 85 

INEUCJ>=0 

GO TO 84 

85 DX=INP(KADJ+1 )-INP(KADJ) 

INEUC J ) -FLOATC INPCKADJ) )+DXX(DS-FLOATCKD) )+ . 5 
MIN = MIN0CMIN, INEUC J) .) 

MAX=MAX0(MAX, INEUC J) ) 

84 CONTINUE 

MAX=MIN0CMAX, MAXAMP) 

C FILL IN TOP BLANKS 

IF ( MAX . EQ , MAXAMP) GO TO 46 
DO 45 I = MAX, MAX AMI 
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45 CALL PLOTITCIPLOT, IPLOT, 13, IZEF'O) 

46 DO 41 I=MIN,MAX 
ICNT=3 

IR = MAX + mN-I 
C FORCE EHD=0 

DO 43 K = 1 , IHIGHLl 

IF((IhfEy(K;-IR):t:(INEU(K.+l)-IR) ,GT.0) GO TO 43 
C PULSE POINT FOUND 
ICNT=ICNT+1 
IPLOT( IpNT)=K-l 
43 CONTINUE 

CALL PLOTITC IPLOT, IPLOT, ICNT, IZERO) 

41 CONTINUE 

IF(MIN,EQ,0) GO TO 70 
C FILL IN LOWER BLANKS 
DO 47 1=1, MIN 

47 CALL PLOTITC IPLOT, IPLOT, 13, IZEF’O) 

70 CONTINUE 

78 LASTC1)=0 

PM = 0. 

FMAX = 1 ./(S..tDXR5tCALTLJD) 

UR I TEC 3, S3) FMAX 

93 • FORMAT (31HLOG10 SPECTRUM . FMAX F9 . E , 6H CVC/M ) 

URITEC3,2S) 

DO 9S I=2,NFFT21 
96 PM=AMAX1 CPM, PAVC I ) ) 

CALL SETQRDCLU, IGRD, 16, ISPAC) 

DO 99 I=1,NFFT21 

IPLOTC 1 ) =1399 . XC ALOG10C AMAXi ( 1 . E-5 , F AV C I ) /PM ) )/5. +1 . )+. 5 
DO 99 J=1,IRPT 

99 CALL PLOTITC IPLOT, LAST, lONE, lONE) 

DELTX=DXRXALTUD 
URITEC3,73) DELTX 

73 FORMAT (30H1LOG OF PROB ( SPACING ), DELTA X= ,F9.5) 

URITEC3,22) 

CALL 3ETGRDCLU, IGRD, 16, ISPAC) 

SUM=0. 

DO 82 1=1,256 
82 SUM=SUM+PDFCI) 

DO S3 1=1,256 
PDFC I ) =PDF( I )/5LIM 

IDAT = 2S0 , X( ALOG10C AMAXl C , 00001 , PDF( I ) ) )+5 . ) + . 5 
DO S3 J=1 , IRPT 

: 83 CALL PLOTITC IDAT, IDAT, TONE, IMl ) 

STOP 

' END 

.U EXT FUNC 

1 IPARM INT4 VAR 

IHEDCK INT4 VAR 

IHED INT4 VAR 

INN INT2 VAR 

■ INF INT2 VAR 
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JAT 

REAL 

VAR 

ISTAT 

INTa 

VAR 

IPLOT 

I NT a 

VAR 

LAST 

INTa 

VAR 

PULSE 

REAL 

VAR 

VAR 

REAL 

VAR 

PDF 

REAL 

VAR 

IGRD 

INTa 

VAR 

INEU 

INTE 

VAR 

TBL 

REAL 

VAR 

Z 

CMPX 

VAR 

CMPLX 

EXT FIJNC 

PAV 

REAL 

VAR 

lONE 

INTa 

VAR 

I2ER0 

INTa 

VAR 

IS 

iNta 

VAR 

ITUO 

INTE 

VAR 

13 

INTE 

VAR 

IE56 

INTE 

VAR 

ISPAC 

INTE 

VAR 

IRPT 

INTE 

VAR 

MAXPLS 

INTE 

VAR 

IBIG 

INTa 

VAR 

I AMP 

INTE 

VAR 

DELTAT 

REAL 

VAR 

NFFT 

INTE 

VAR 

-U 

INTE 

VAR 

YES 

REAL 

VAR 

DTOR 

REAL 

VAR 

DLT300 

REAL 

VAR 

16 

iNTa 

VAR 

I M 1 

INT2 

VAR 

SYS 10 

EXT 

FUNC 

FOUR IT 

EXT 

FUNC 

NFFTEl 

INTa 

VAR 

40 

LABEL 

@H 

EXT 

FIJNC 

50 

LABEL 

PRPULS 

REAL 

VAR 

4 

LABEL 

5 

LABEL 

30 

LABEL 

31 

LABEL 

ANGLE 

REAL 

VAR 

BEAM 

REAL 

VAR 

ALTUD 

REAL 

VAR 

ANGRAD 

REAL 

VAR 

BMRAD 

REAL 

VAR 

BM2 

REAL 

VAR 

ANGMIN 

REAL 

VAR 

A MAXI 

EXT 

FLING 

ANGMAX 

REAL 

VAR 
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JXR 

REAL VAR 

SIN 

EXT FUNC 

COS 

EXT FUNC 

DANG 

REAL VAR 

42 

LABEL 

I 

INTE VAR 

GO 

LABEL 

1 

LABEL 

2 

LABEL 

IFILE 

IHT2 VAR 

7 

LABEL 

IBLKS 

INT2 VAR 

IRECS 

INT2 VAR 

IDOCS 

IMT2 VAR 

IDATS 

IHTE VAR 

LO 

INT2 VAR 

IHI 

INT2 VAR 

I DOCS 1 

I NT 2 VAR 

IREC2 

INT2 VAR 

IDATS2 

IMTE VAR 

MAX0 

EXT FUNC 

NAX02 

EXT FUNC 

NINO 

EXT FUNC 

MIN02 

EXT FUNC 

IHIGH 

I NTS VAR 

IHIGH2 

I NTS VAR 

( IHIGHLI 

INT2 VAR 

SAMPS 

REAL VAR 

. U 

EXT FUNC 

DUR 

REAL VAR 

IBEG 

I NTS VAR 

KTOT 

INT2 VAR 

9 

LABEL 

- nil 

I NTS VAR 

128 

LABEL 

80 

LABEL 

H as 

LABEL 

228 

LABEL 

81 

LABEL 

BVTOHU 

EXT FUNC 

100 

LABEL 

17 

LABEL 

1;. K. .. . 

I NTS VAR 

|; .X; 

. REAL VAR 

j ; 200 

LABEL 

i; SUM 

REAL VAR 

!:' 8 

LABEL 

1 XMAX 

REAL VAR 

; 10 

LABEL 

i • 6 : 

LABEL 

i 3 

LABEL 

i ao 

LABEL 




ai 

23 

SETGRD 

12 

lEDAT 

i^LOGie 

i, V 

U 

PLOT IT 
23 

N 

^TOT 
;si 
i. R 

IKSTART 

'DELTAR 

iDELTAS 

Ikend 

IIBI 

|PD0 

!51 

K2ER0 

iS2 

iXZERO 

^64 

I3CALE 

;S2 

ISQRT 

'72 

DRSTRT 
;K0 
i IFIX 
|DT0T 
'FLOAT 
: FL0AT2 
:25 
i KSET 
‘ 4S 
I XR 
: RNS 
' DR 
: Ds 
! KD 
! KAD J 
i 58 
I 49 
1 57' 

I ANGTOT 
I 6 6 
i 67 
i ANG 


LABEL 
LABEL 
EXT FUNC 
LABEL 
INTE VAR 
EXT FUNC 
EXT FUNC 
INT2 VAR 
EXT FUNC 
LABEL 
REAL VAR 
REAL VAR 
LABEL 
EXT FUNC 
INT2 VAR 
REAL VAR 
REAL VAR 
INT2 VAR 
LABEL 
REAL VAR 
LABEL 
INT2 VAR 
LABEL 
REAL VAR 
LABEL 
REAL VAR 
LABEL 
EXT FUND 
LABEL 
REAL VAR 
INT2 VAR 
EXT FUNC 
REAL VAR 
EXT FUNC 
EXT FUNC 
LABEL 
INT2 VAR 
LABEL 
REAL VAR 
REAL VAR 
REAL VAR 
REAL VAR 
INTE VAR 
I NTS VAR 
LABEL 
LABEL 
LABEL 
REAL VAR 
LABEL 
LABEL 
REAL VAR 
LABEL 


b‘9 

LABEL 

79 

LABEL 

65 

LABEL 

I.H 

EXT FUNC 

89 

LABEL 

77 

LABEL 

MAX AMP 

INT2 VAR 

MAX AMI 

INTE VAR 

TIMING 

REAL VAR 

ITMINC 

INT2 VAR 

70 

LABEL 

7S 

LABEL 

MAX 

INT2 VAR 

MIN 

INTE VAR 

71 

LABEL 

34 

LABEL 

ALF 

REAL VAR 

ATAN 

EXT FUNC 

BX 

REAL VAR 

ABB 

EXT FUNC 

as 

LABEL 

ANTFAC 

REAL VAR 

67 

LABEL 

RADJ 

REAL VAR 

91 

LABEL 

DX 

REAL VAR 

TOPEAK 

INTE VAR 

97 

LABEL 

95 

LABEL 

IPKN 

INTE VAR 

94 

LABEL 

»P 

EXT FUNC 

98 

LABEL 

CABS 

EXT FUNC 

S4 

LABEL 

85 

LABEL 

46 

LABEL 

45 

LABEL 

41 

LABEL 

ICMT 

INTE VAR 

IP 

INTS VAR 

4i 


47 


Pti 

REAL VAR 

FfWX 

REAL VAR 

93 

LABEL 

96 

LABEL 

99 

la!bel 

DELTX 

REAL VAR 

73 

LABEL 

82 

LABEL 

83 

LABEL 


/ 



; 5 EXT FUNC 

I ,V EXT FUNC, 
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Sample Printout 


Program ECKPDF 


Any^litude Probab^l vtj^es 


^ ‘ A •: * 
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Program ECKPDF - Sample Printout 
Average Time Pulse 


AMPLITUDE 



RANGE 


i'OLDOUT niAMB 



Program ECKPDF - Sample Printout 
Average Pulse in Range 
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rOLDOUT FRA^^rH 






r 





h / 


III .in I rill 1.1 III ' "u. u v »4 i,i; 11 , POWER 




FREQUENCY — ♦ 





1 



Program ECKPDF - Sample Printout 
Average Power Spectrum 



POOR quality fX)LDOUT FRAivfH 




>.1^ 




Linear Time 
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L|,near Range 




Linear Angle 
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Appendix F 

SHORT PULSE RADAR SIMULATION 
Program Listing and Sample Printouts 


Program SIMUL - Program Listing 


PH«iE 0001 




tflSSM 

SIl-WL PROG a D SHORT PULSE SLmTTER SIMULATION LDD Z/7(. 

EPLST 

C PROGRAM TO SIMULATE A M-41'IFICAT ION OF EON ( 4 ) IN THE TIMF. 

C DOMAIN OF -SPECTRUM OF POWER “ LEVINE OF AUG 1974 

C THE STATIONARt POINTS ARE FOUND FOR A GIVEN SURFACE 
C USING FFT TECMNICiUES 

C BOTH A SINUSOIDAL AND A GAUSSIAN SURFACE 
C COMPONENT ARE ALLOVitD. THE SURFACE 
C IS GENERATED IN THE FRFOUENCY DOMAIN 
C TO GET THE DESIRED SPF-'TPUM. THE 
C RANDOM SURFACE IS PTERt E-MOSkOWIT 2 
C WITH GIVEN PEAK FREQUENCY AND 
C RMS VALUE AS INPUT. THE SINUSOIDhL 
C PART HAS THE SAME PERIOD BUT ANY 
C OTHER CHOSEN RMMS VALUE. 

C PROGRAM BY TECHNO-SCIENCES. INC 137C 

C 

C 

If'Y^LICIT INTEGER*^ (I-N) 

INTEGER*4 IRAN 

COMPLEX 2( lOr.4 ZD< 1CE4 n 2DD( 1024 CFF-LX. CON JG 
COMPLEX U. 22 

DIMENSION TBl< S13 )- IGPD' 21 ). IPLOTt 21 )• LAST'. 31 ) 
DIMENSION RANGE< 1024 '.CROOK 1U24 ) 

EQUIVALENCE < RA'NGEC 1 2^ 1 ) J. ( CROOK 1 }• HDD( 1 ) ) 

TAN( < J»SIN< X VC OSt X ) 

C 5PEH X )-PIERSON-MOSk0UIT2 SPECTRUM 

SFEC( X )-A2vEXP( -AMINK 130. . eC2U4/< X*X > ) )/X ) 

C X-ANGULAR SURFACE FREQUENCY 
C I RAHa FANDOM NUMBER '‘aENERATOR 
C GAUSSIAN PULSE SHAPE 

PULS£( X >EXF( -( t. X.'SIGFLS )-»*2 l-'Z. ) 

C ^ SIN(. X )/X )»*2-ANTEMNA POWER'iAIN PATTERN 

GAIN< X SIN( X/'ANGVAL )/AMAXl< ASS( X/ANG^AL ). 1. E-3) )**V2 
DATA SIGTIM.". 00*/ 

DATA IRAN/ 1234321.x 
DATA NFFT/ 1024.' 

DATA LU/C5." 

DATA P2/'5. 2-S31S53/ 

DATA C. DT. PULSTM. FO/300. . . 00:25. .01. 13900. X 
C C-SPEED OF LIGHT. DT-TIME INCR. PULSTM-PULSE DUPAT10N 
DATA 13/3/ 

DATA 14/4. 

C DT0R -DE'iREES TO RADIANS CONVERSION 
DATA DT0R, . O 17*5323.'' 

DATA lONE. IBIG/'l. 30000/ 

DATA IGDPTS/7/ 

W‘O'-F0*P2 
NFFT 2» NFFT /'2 
NFFT 21 -NFFT 2* 1 


■ QVALrri 
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WPITE^O.25) 

2S P0Rr'IAT< 59HENTEP RI1S 5INU30IDAUM. ). 

IRMS RhNDk)^X P-M->f M. J.PERIfoDiM. KFA) ) 

READ(0. 24; RMSSIN. RM3RAN, PERIOD 
24 F0RMAT<8F3; 

WRITE( 0. 27 ) 

27 FORMAT ( 33HENTER RC 3DB CUT0FF*PIJLSE DUR(F3) ) 
READ^O.24; FOTS 
WRITER 0.31 ; 

31 FORMAT < 30MWAVE REPEATS. PULSE REPEATSdS); 

PEAIUO. 32; IRPTW. IRPTP 

32 FORMAT <315; 

WRITER 0. 4 1 ) 

4 I F0RMAT( 1 IHBLOCF 37< 15 ; > 

READC 0. 32 ; I6LK3 

PH0-EXP< -P2*F0TS*DT 'I’ULSTM; 

RH01-SOKT< l.-RH0#RHO) 

BG2U4-P2*Pc* 1. S/T PERI0D*PERI0D ) 

A£«P2*RM5RAN*RMSRAN*BG2U4 

ALTUD-2000. 

C2-C-^. 

C CON-ERT FROM PUl.SE TIT« SIGMA 10 RANGE SIGMA 
SIGPLS-3IGTIM#C2 

C RNGINC-ASSUMED PULSE START. END POINT 
FNGINC-4. *SIGPLS 
ANGLIN- 18. 

ANGCTR- ANGLIN* DT OR 
BMWITH-24. 

ANG3DS-BMWITH*DT0R/2. 

ANGMAX- ANGCT F f ANG3D6 
ANGMI N- ANGC T F - ANG3D8 
ANGVAL»ANG3DB.'l. 33 
lfO-ALTUD*TAN( ANGMIN ) 

C SET UP GRID 

GDINC- 1333. .'Fl.OATC IGDPT5- 1 
D0 5 I- 1. IGDFTS 

5 IGRD< I ;-FL0AT*. I- 1 ;*GDINC +.5 
t3DB-ALTUD*< ThN( ANGMAX ;-TAN( ANGMIN ) ) 

Dir-t3DB. FLOAT', NFFT ) 

DyZ-Dt/2. 

•<02-T0-DY2 
C DF*Dt-l.NFFT 

DF- 1. /Y3DB 
RTDF-SQRT( DF . 

DW-P2*DF 

C DW-DELTH OMEGA. DF-DELTA FREQ 
C SET UF FFT TABLE 

CALL FOUR 1T<E. NFFT. TBL-0) 

C SET BLOCK COUNTER T0 2ER0 BEFORE LOOPING 
ICOUNT-0 

C SET UP RANDOM FART 0F WAVE TRANSFORM WITH SPECTRUM SPEC< X ) 
C DC PARTS-0 


i 
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z< I )-*CMPi.x>: 0. . Cl. ) 

2D( 1 >-CtiPLX( 0. . 0. ) 

HDD< 1 )»CMPLX( 0 . . 0. ) 

D0 1 I-2.NFFT21 
X-FL0AH I- 1 >*DV/ 

I J-CMPLX< S&RTC FFECC X ) IRftN ). 0. ;*RTDF 

1 C0NTINUE 
IFREQ-T3DB/PERI0D*. 5 

C SET UP M0N-RftND0M V/FIVE C0MP0NENT 

2t IFREO )-2< IFREQ )fCrFLX< RMSSIN.^QRTt 2. ). 0. ) 

C SET UP F0R PLidT SCALE FhCT 0RS 
3CALE»RM3SIN**2 
X«FL0AT( IFREt- 1 )*DW 
SCHLE1-SCALE»X4X 
SCALE2»SCALE1*X»X 
D0 2 I-2-NFFT21 
X-FL0ATI 1- 1 J*PVr 
XX-SPEC( X )*£. *PF 
SCALE- SCALE* XX 
XX-XX*X*X 
SCALEl-SCALEl+XX 
XX-XX»X*X 
SCALE2-SCALE2+XX 
W-CtiPLX( Q. . X ) 

2DC I >L*2( I ) 

2I'D( I )»W*2D( I ) 

2 CONTINUE 
SCALE-SQRTt SCALE ) 

SCALE 1-SQRT< SCALE 1 ) 

SCALE2-SQRTC SCALE2 ) 

D0 3 I-2.NFFT21 
IC-NFFT+2-I 

2< IC )-C0NJGt EC I ) ) 

EDC IC )-C0NJG( EIK I ) ) 

2DDC IC )-C0NJGC cDDC I ) ) 

3 C0T4T INLIE 

CALL F0URITC E> NFFT. TBL. I0NE ) 

CALL F0UR1TC ED. NFFT. TBL. I0NE ) 

CALL F0UR1TC EDD. NFFT. TBL. I0NE ) 

C FIND STATI0NAFY F0INT3 ASSUMING LINEAR INTEFP0LATI0N 

C AND PL0T E(t J. FIRST e.2D DERIVATIVES 

CALL SETGRDC LU. IGRD- IGDFTS. IBIG > 

TMDIF-C T02+DX ).'ALTUD - REALC HDC I ; ) 

K-0 

D0 4 I- 1. NFFT 
X-FL0ATC I ,^*DT + 102 
TANTHT-t.-IALTUD 
TANALF-REALC ?D( I ) ) 

TANDIF-TANTHT-TANALF 

IF C TANDIF*TMDIF. GT. Q. ) G0 TO 6 

DLT-TMDIF/C TMDIF-TANDIF ) 

DLT 1- 1. -DLT 


ORIGINAJUPAr . 
OF, POOR QU. ' ' 
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r«t*DLT*DY 

HEKiHT-REHLC E< 1^1 J )*DLT l + DLT*h!ErtL( 2t I ) ) 

DERIV-RErtL( 2Dt I- I ) )*DLT l + DLT*THNALF 
riERIV2«REAL( EDI'< I- I ) i*DLT l*DLT*REflL( 2DDf I ) ) 

K»K* I 

R^N(iE^ K )-SQRT< 1*Y*( rtLTUD-HEIGHT )**2 ) 

MN^i-ATHN^ DERIV ) 

KHDCRV-( ( 1. ♦rERIV«DERI¥ 1. 5 J^PERIVE 

CR00T( K )-3CRT( ftBSt RADCRVA < RhDCRV-RHNGE^ K > ;#RANGE» K ) ) ^ ) 
CR00T( K )-CR00U K )*GRIN( rtNG-ANGCTR ) 

6 C0NTINUE 

IF( IRPTW. EQ. 0 )G0 T0 4 

IPL0T( 1 )-20. ♦REPIL( H( I ) >/SCPLE + 233. 5 

IPL0T(2)-2O. tREftHaDd J )/SCALE !♦ 700. 

IPL0T< 4 )»20. ♦TftNTHT^'SCrtLEl + 700, 

IPL0T( 3 )-20. ♦REAU 2DD( I ) J/SCHLE2+ U67.5 
D0 30 J-l.IRPTV/ 

C REPEAT T0 SPREAD SCALE 

30 CALL PL0TIT( IPL0T. LAST. 14. I0NE ) 

4 TMDIF-TANDIP 

NUM-K 

C REORDER IN INCREASING RANGE 

10 D0 7 K-2. NUh 

IF< RANGEC K- I ). GT. RANGE! K ) ) G0 T0 H 

7 CONTINUE 
G0 T0 S 

3 X-RANGE(K) 

RANGE! K )-RANGE( K- 1 > 

RANGE! k- 1 )«X 
X-CR00T! K ) 

CR00T! K .I-CR0CT! K- I ) 

CR00T! K- 1 >X 
G0 T0 10 
3 C-’'NTINUE 

c generate outputs now 

TIME-! RANGE! 1 ;-RNGINC VC2 

Y-0. 

LAST ! 1 )-0 
IPL0T! 1 )-0 
LRITE! 3.21) 

21 FORMAT! IHO) 

11 TIME-TIME^DT 
RG-C2*TIM£ 

C RG-RANGE FOR PULSE CENTER 
C RGMAX. RGMIN-STAPT END 
RGMAX-RG+RNGINC 
. RGMIN-RG-RNGINC 

C FIND REFLECTOR LIM1TS.MINV.MAXV 
C I. E. THOSE IN THE RANGE BIN 
D0 12 I-l.NUM 

IF! RANGE! I ). GE. RGMIN ) G0 T0 13 

12 CONTINUE 


& 


Original page p 

OB Poor qualiV.' 


13 

IC0UNT-IC0UNT* 1 

IF< IC0UNT. LT. lELKS ). G0 T0 40 

3T0P 

MINV-I 

14 

D0 14 I-MINV.NUM 

IF( RANGE( I ). GT. RGMAX ) G0 T0 15 

C0NTINUE 

15 

I«NUM-» 1 
MAXV-I- 1 

2Z-Ct>lPLX< 0. . 0. J 
C GENERATE 0UTPUT3 

C 6. PL0T THEM 

lo 

IF( MAXV. LT. MINV )G0 T0 20 

D0 16 I-MINV.MftXV 

ANG-W0*( RANGEC 1 )-RANGE( MINV ) )/C2 

W-CMPLXC C0S( hNG ). 3IN( ANG t ) 

22-22*CR00T< I ;*PUL3E( RANGEC I )-RG 

C0NTINUE 

C FILTER 

20 

X-RH0*t*RH01*t. CA8SC 22 )**2 ) 

24 

t-X 

IFC IRPTP. EG. 0 JG0 T0 11 
IPL0TC 1 )• 1400. ♦XVALTUD 
D0 24 I- 1. IRPTP 

CALL PL0TITC IFL 0T. LAST. I0NE, I0NE 

18 

F0RMATC IH .4E 16. 6. 417) 

. U 

IRAN 

G0 T0 11 
END 

EXT FUNC 

INT4 VAR 


CMPX VAR 

2D 

CMPX VAR 

HDD 

CMPX VAR 

CI^LX 

EXT FUNC 

C0NJG 

EXT FUNC 

U 

CMPX VAR 

22 

CMPX VAR 

TBL 

REAL VAR 

IGRD 

INT2 VAR 

IPL0T 

INT2 VAR 

LAST 

INT2 VAR 

RANGE 

REAL VAR 

CR00T 

REAL VAR 

TAN 

STATE FN 

. Q 

EXT FUNC 

. P 

EXT FUNC 

XB 

F0RM PAR 

SIN 

EXT FUNC 

C0S 

EXT FUNC 

SPEC 

STATE FN 

xa 

F0RM PAR 

A2 

REAL VAR 


Driginal page 

V?J POOR QUALI 


PA(JE OOOto 


EXP EXT FUNC 

rtMINl EXT FUNC 
BG£;U4 REAL VAR 
. R EXT FUNC 

PULSE STATE FN 
XI F®RM PAR 

SI'iPLS REAL VAR 
GAIN STATE FN 
XH F0RM PAR 

ANGVAL REAL VAR 
AMAXl EXT FUNC 
ABS EXT FUNC 

SIGTIM REAL VAR 
NFFT INT2 VAR 
LU INT2 VHR 

F2 REAL VAR 

C REAL VAR 

DT REAL VAR 

PULSTM REAL VAR 
FO REAL VAR 

13 INT2 VAR 

14 INT2 VAR 
DTBR REAL VAR 
IflNE INT2 VAR 
IBIG INT2 VAR 
IGDPrS INT2 VAR 

REAL VAR 
NFFT2 XNT2 VAR 
NFFT21 INT2 VAR 
25 LABEL 

C-H EXT FUNC 

2o LABEL 

RMSSIN REAL VAR 
RMSRAN REAL VAR 
PERI0D REAL VAR 
27 LABEL 

F0T3 REAL VAR 

31 LABEL 

32 LABEL 

IRPTV/ INT2 VAR 
IRPTP INT2 VAR 
41 LABEL 

lELKS INT2 VAR 
RH0 REAL VAR 

RH01 REAL VAR 
SORT EXT FUNC 
ALTUD REAL VAR 
C2 REAL VAR 

RNGINC REAL VAR 
ANGLIN REAL VAR 
ANGCTR REAL VAR 
BMWITH REAL VAR 


ORIGINAL PAGE IS 
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rtNG3DB REAL VAR 
ANGMhX real VAR 
ANGMIN REAL VAR 
to REAL VAR 

GDINC REAL VAR 
FLOAT EXT FUNC 
FLOAT 2 EXT FUNC 

5 LABEL 

I INT2 VAR 

. T EXT FUNC 

t3I'B REAL VAR 
DY REAL VAR 

DY2 REAL VAR 

Y02 REAL VAR 

DF REAL VAR 

RTDF REAL VAR 
DW REAL VAR 

FOUR IT EXT FUNC 
IC0UNT INT2 VAR 
40 LABEL 

»P EXT FUNC 

1 LABEL 

X REAL VAR 

RANDG EXT FUNC 
$M EXT FUNC 

IFREQ INT2 VAR 
tL EXT FUNC 

»A EXT FUNC 

SCALE REAL VAR 
SCALE 1 REAL VAR 
SCALE2 REAL VAR 

2 LABEL 

XX PEAL VAR 

IK EXT FUNC 

3 LABEL 

IC INT2 VAR 

SETGRD EXT FUNC 
TMDIF REAL VAR 
REAL EXT FUNC 
K INT2 VAR 

4 LABEL 

Y PEAL VAR 

TANTHT REAL VAR 
TANALF REAL VAR 
TANDIF REAL VAR 

6 LABEL 

DLT REAL VAR 

DLT 1 PEAL VAR 
HEIGHT REAL VAR 
DERIV REAL VAR 
DERIV2 REAL VAR 
ANG REAL VAR 


PAGE OOOS 


AT AN 

EXT PUNC 

RADCRV 

REAL VAR 

. A 

EXT FUNC 

30 

LABEL 

J 

INT2 VAR 

PL0TIT 

EXT FUNC 

NUM 

INT2 VAR 

10 

LABEL 

7 

LABEL 

3 

LABEL 

9 

LABEL 

TIME 

REAL VAR 

21 

LABEL 

1 1 

LABEL 

RG 

REAL VAR 

RGMAX 

REAL VAR 

RGMIN 

REAL VAR 

12 

LABEL 

13 

LABEL 

. S 

EXT FUNC 

MINV 

INT2 VAR 

14 

LABEL 

15 

label 

MAXV 

INT2 VAR 

20 

LABEL 

lb 

LABEL 

CABS 

EXT FUNC 

24 

LABEL 

13 

LABEL 

. V 

EXT FUNC 

0000 ERR0RS 


QUALTiyj 


Program SIMUL - Sample Printout 


79 


Simulated Surface Second Derivative 


Intersection Defines Specular Points* 



Simulated Surface First Derivative 


Simulated Surface 
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Appendix G 
SUBROUTINES 
Program Listings 


Tit FU«&1 


program PLOTIT 




84 


PPO'l* PLOTI’ 


C j-OciFJ'.'riii 


. o^nooi 
‘^OOOC'L'I 

■jC'C'OuOI 
•C CnJC 41 
i-C O X Z 1 
vC IGuCX 
‘ 'CiO'.OI 
C'-' jO'.cI 
0 i.G 0 ; 4 f 
;0 : 


I'tPU 

AiF.F 

•:afo 

SC >0 
z?F.ei 
IZZ' 

4ir0 40C0 PCOCF 


JI‘S 

QCCj 

fffr 

ii;o 


LL.-OZoi 

E 170 
24E0 


u.'OC^dl 

PSCl 


■ UC024I 

CSDO 

OCAC 

'■•roocSi 

SCcD 

i210 

■ - -Ti'V-CJ 

SOCD 

315 7 

coooBoi 



: COO/ 21 

|21S 


'X007.4I 

Dicr 

OCC 4 

Luoor SI 

4iiCC 

OC C .' 

OOuOCCI 

I'2u0 

•ilCS 

0000401 

D2C0 

.31 IS 
OCC 5 

0000441 

4.2EE 

00 

1 lEl 


'.0004mI 

2TE2 


0'.-- 

21 IE 


Oc004LX 

43CD 

4EOO 

0C-30S4I 

2312 


C'.'00S'4hI 

OOOobc.1 

c4:d 

C5C0 

GS;3 

OOOOSCI 

< 1S3 


OuOOSEl 

gcoocii 

<: *on 

0577 

TSCO 

1 L e 


• O'.'OUe^l 
••OGO'Jal 
'.'■■OOta 


ijOOC VI 
OOOO'EI 
00i:i0oii 
nOOOciel 

IjvOO^lil 

OCOOKI 

uoooaci 

OOOOSJI 

OOOG&CI 
OOOO.-EI 
1 -' OOiJi'iil 
'OOOH 4 i 
ijOjOt-ai 
POO'JhCI 
COOr'BDI 
CC'uhtl 
• '■ jijfiil 
COCOShI 
OuouL';i 
'.'OijOCOI 


E20E 

4£FF 

SOFO 

i 7 Fl 

EOFO 

4300 


I’OrtO 

4&EF 

■:4F0 

50F0 

i 7 En 

iil 9 

46 FF 
4<:E£ 
I lEi 

sseo 

Z4B I 
ECbO 
ESBO 
42T0 
24 B 0 
5C60 

41 ;o 

t’ti. 

42 10 
46BC 


OCCO 
SI 14 


21 •£ 
iC'GO 


iiDEE 

OCCO 

FFFO 

ACf 

FFM 

OCC 4 

OOCO 

OC'CO 

dCE 4 

eoEo 

3C>r-2 

aciE 

dCIc 

302c 

4ECJ 0000 


3 

-.0 


.2 

;4 


.-4 

rS 


r. 

"0 

i 


73 

72 

40 


■ -3 

<13 

'0 


?CP>»T 

ENTPIE3 AT 3ET<}PI’ AHP FLe>riT 

^0pi4an callable, h't '.'E called 

FIRST *0 SET UP 3FI7 FlITTIN';. ChLl 13: 

CALL SEI'.PI". ..O. ICPID' 1 . HGR ID. M -SPACE ) 

WHERE I 


L0*PL'3TTEP L. 0. 

I<iRID li ‘ 

M'lPID.tHE NOI^EF iF THEiE POINT : 


AN AFRAlt CP e iPIP c3CHriAflS 


ti.<:.PHCE.*NijfieEp OF pt-n.TS fft./eeh -iPin 
LINES IN ’ME .-DIFECTIiS. IF -0. MO IPIC lINES 
T0 plot, call 

CALL PLCTITi lOATP. .MSTrviH./lPaiNTS.NCONNECT ) 
WHERE UATft 13 H CATS mPPHT 
NCONNECT. 0 TO FlVT POINTS CNLt. 

NZ T0 ccnne:t T-fiE ffi«fnt point with 
THE LAST. IF .L', -:. ■•NMEC* FR*)H Oi 3flP TltPf i. 
NPflINTi.NUf'EEP JF POINT: ill DATA APPB«T 
LASTDATA is an AFPAT *0 HOI D ’HE -AST 
POINT; OF Dir-€NS:CN SOUAL TO MPOINTS. 

LASIDATA 13 ijSFD EVCM IF ;in C-fUNECTIBN 
IS I'EsIFED. AND I1JST PE iI»E’.. 

ALL values -PE «-PY’7 rUTEOERt. 

AND APE ASSUTED T* _IE BETWEEN 0 W 
MA. FT 3- 1.133?. Thcy apf LIMITED IC 
this PAiiGE if THET ape no* on ENTPV. 


SE' 'iPD 


iPP 


IP’S 


ENTFT 

EXTRA 

STM 

LH 

NHI 

ST 

SIS 

ez 3 

LHI 

eal 
DC. 3 

til 


SETCPD.PLaTIT 

. 0 

10. FEB 

14. L< l5 i 

15. l-'-PFC 
IS. P 15 

14. :o 
CCS 

1 1 . c ■ ■ 

15 . . 0 

E170 
14.0 
: c. 


Sh'-E peg 3 

■iPT Nure ;;f ap.,s. 

i1a;k is 70 F. 

- SAVE 
NOW .TP APGS 


ScMD EFP TE35 
■- -iUIT 

SET JP GRID APSAt 



PI 

ST 

14.G0'JI' 13 ■ 

EEP0 This one 

.14 


ST 

12,GPID< 13 • 

HI.L *NE5 F0F ^EP’ICAI. LINE 



SIS 

.3. 4 

DECP CNT 

4© 


ENDIS 

PI 


47 


LIT 

12. 4I IS ■. 

ir7T ARCS 

40 


LA 

12. C( 12 > 

SET lU 

•*3 


3T8 

12, lPITE* 1 


50 


STS 

12.LPGRID.1 




LM 

14, C< 14 

GET NUTIEEP ‘JF grid LT.NES 



SLLS 

14. 1 

•2 F0R ACCESS 

t j 

P* 

SIS 

14. 2 

DECR PTP 

§4 


Bfi 

P3 

OUIT WHEN - 

■ n 


LH 

12. 0< 13. 14 t 

I.ET NEXT F0IMT 



eNHis 

P4 

OF LlfllTS 

cr • 


LIS 

12.0 

SET.Q 

t. . 

P4 

CLAI 

12.MA..PTS 

■;f T00 la»ge 

59 


BLS 

P5 


•40 


LMI 

12. MA..PT-:.- 1 

limit UPPER 

1 

P = 

SBT 

12. GOUT 

SET GRID P0INT 

<C2 


B 

P2 

.E I NE ‘ T 

6Z 

P3 

LH 

15. IS- 

jcT -t SPACES 

t'4 


ST 

IS. NXSPACE 


tb 


SIS 

15. 1 

f 2RCE HXIS HT JEGifiN .NG 

« 


ST 

15. -TCCUHT 


r7 


8 

P 13 

E TIT 


• plotting ehtfy point 


;i PLQJTIT 

stm 

IOpFL^ 


■ c 

LH 

14. 0M3 J 

iET hRGS 


NHI 

!3 p r^FFC 

Fix UP PIS 

: 4 

ST 

;s . F :e 


; c, 

SIS 

14. 1 j 

Ck NUrtE hPG; 


SNZ 

EP® 

:r N0N SEP, EM' CN 

‘7 

LM 

U. 4f . 5 ' 

GET afg list 

'S 

LH 

15. 0‘ 15 ' 

CONNECT <ALUE 

73 

LH 

14. 0< 14 j 

NLltiER 0F PCIN’S 

■jO 

SLLS 

14. : 

• 2 F0R ENTRY 

61 

L 

: 1. -CJ-.NT 

Ck f-3R A GPI. line 

•i2 

AI3 

i 1. 1 


l :3 

ST 

I l. :iCJ‘JNT 


•.*.4 

CL 

I 1. 


■:-b 

ENE 

P14 



lIS 

i 1. 

RESET COUNT 


ST 




SvC 

1 . r* 

..rite -EFTI-.-.l lit. 

^16 

SIS 

14. £ 

NOW -ABVE T0 l.m 3T r 

•■-0 

EM 

PI.: 

until D 0NE 

31 

lH 

l 1. C' l4. 14 , 



'ESS 


ORIGINAL' PAGP iq 
POOR Qt® 
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tumxim 


ROQRhM TI6 PLAT 6«Tft Lt'D l.rc 



CiOdoEc 
ocooee: 

OOOOF? 
-•000F4 
jijOl IF A 

ooooFfi: 
oocOFc: 
OOOOFE, 
•jijo 1C4: 

000 IfiMi 

oooicic: 

000 1 oe: 
ocoi 12: 

0001 ;j: 
0001 '..i: 

‘•001 !c: 

IJOOIOO. 

onoi;.£: 

000 lr4: 
000 i4s: 
‘■'00 1 
000 i2c: 
•oooioe: 
000 i:o: 
'jcoi;:2: 
■-'OOIO';': 
0'?013m1 
OOOlOEI 



4EC0 0000 

OChC 
SMi 
ii. 12 


AOf 2 

4EC0 0000 


0S7A 

0E>7 


4E00 0000 
4EC0 OCQO 


0^/S 


057- 

tECJ 0000 


4g8D 

>.£80 
5£HP 
50h8 
2 764 
22 IS 
27E2 
4210 
4tBC 
2J12 
2460 

'iseo 
21 A3 

■:eeo 

OeFF 

Siief 

4A8I' 

2312 

2480 

'.5h0 

2 1 A3 
0380 
4C‘8I> 
0568 
2'3w4 
0€F8 
03HP 
026F 
■'SmO 
2»h 1 
05BH 

22 -4 
4300 
e I n 
D180 
43CF 


0000 QCFF 
0000 '3573 

210U 

0000 

0000 '32ECI 
0000 03581 
OCQO OOCO 
0000 OOCO 
2100 
0000 

0000 013CI 
0000 02381 

0000 :cco 

0000 OOCO 


31C0 


FFF3 

4 CC 8 

3C2E 

0C14 




STH 

"iZ 


6 

i'4 

P14 

LHI 

55 

PIS 

L 

' s 


ST 

■?7 


SIS 

->a 


ENMS 

.4 

P8 

*^1 3 

U'U 


6M 



LH 

zr,l 


BNl-5 

11'3 


LIS 


104 

105 
lOb 
;c>7 
loe 

.03 


P8 


PS 


i.LMI 

6LS 

LHI 

LP 

6P5 

eris 


1 .0 


STH 

1 ; 1 

CONNECT 

-H 

! 2 


bNMS 

1 : 

BAP / 

uI5 

1 i4 

PIO 

CLHI 

; ;5 


SL3 

i . ^ 


LHI 

\ ^ 

PI 1 

STH 

1 • ,» 


CLP 

'• -1 


ENLi 

Ir.O 


LP 

1..1 


LP 

U2 


UP 

123 

P12 

5BT 

12.» 


AI3 

125 


CLP 

tiS 


ENLS 

:er 


6 

1 c 5 

P7 

s.-: 

.2 3 

P13 

LPI 

1:0 


6 

« - t 

« 



LINE5IHE 

eou 

173 

Mfi.:FT3 

eoy 

’ :*i 


ALIOT 

:t5 

WRITE 

DCX 

1 **« 


OAC 

137 

WRORID 

DCX 


133 




r(tr,r. 

1 ; . 0^ 1 3. 1 4 > 

PU 

ll.LlliE5r2F.-3 

10. ( 1 i I 

IC.liUI' 11 

11-4 
F 15 
14.2 
P-' 

I i.O' 12. 14 > 

PS 
1 1.0 

1 1. rWXPTS 

1 1. Mfl ;?T3- 1 
15. 15 
CONNE'l r 
EmP 

1 1.0< 13. 14 > 
lO.C" 13. ; 4 J 
PIO 
10.0 

10, MA p’5 
PI 1 

;o. 1 

1 1.0( 1 :. 14 ' 

1 1 . 10 
P12 
15. 10 
10 . 1 1 

11. IS 
10. OUT 
10. 1 
11. 10 
P12 

P4 

1.V.PITE 
10. PE". 

20( 15 , 

l.'S 

1400 

4 

2100. .0 


OUT.'5UT*LIMESI2E- 1.0.0 

2100.0 

'IRID. '-.?ID*LI(4E3I5E- 1-0.0 


SET OF 61 J I PUT lUFF 
uITh OPID UTNFS 

OECR 'ZOOUT 

iv?U DO DhTA P5HAT5 
UHTIU D«UE 
iET POtNT 
•Zk Llf-UTS 
uIMIT TO 0 TO 


MAKPTS- I 

■Z. FOP C0NHEC1 

NH ItANS CONNE ZT 

IT -.FOPiZE 0 ■- INNEO T' eAp- I TPI. ) 

Fi.it F'iPCE i-IO ZOkT.EC b ( 3T «i>t 

..-T last FOIMT 

■-K UMTS 


STBPE !<EU POINT Ik L»iT 
P .'PCS ' 10 I LT 11 ; 

r-‘'.'»E 


SET SI'S 

uUTIl 

■ 10 '-111' 

'JET NEaT 
liPITE IT OUT 
E-IT 


BtTES 

BITS 


i33 

PEO 

DS 

20 

!40 

PIS 

PS 

4 

141 

NXSPACE 

DS 

4 

142 

X COUNT 

DS 

4 

143 

OPIO 

DS 

LINESIsE* 

144 

UOUT 

DS 

LINESI2£» 

145 

14o 

OUT 

DS 

END 

LINESIsE* 


00 POOR 





Program RANDG 


PHGt 0001 


tfiSSM 

RANDG PRUG GAUSSIAN F:AND0ri NUMBER GENERATOR LDD 3/76 
ERLST 
$FORT 

FUNCTION RANIG(IY) 

C It IS AN INITIALISING VALUE, ODD. LI 9 DIGITS, 

C INTEGER»4. NORNiAL NUMBERS GENERATED BY SUMMING 
C IS UNIFORMS IN I PE USUAu fASnlON. kANDU ROUTINE 
C IS USED (NOT CALLED). 


RANDG 
RANDG 
. Q 
. P 
lY 
X 
I 
I 

FLOAT 


X-0. 

DO 1 I- 1, U 

IY-IY«6553S 

IF( It. GT. 0 ^G6 10 2 

IY-IY+2147483647^ 1 

X-X*FL0AT( lY , 

CONTINUE 

RANDG-X/2 14748364 7. - 

RETURN 

END 

FUNC/^UB 
FUNC VAR 
EXT FUNC 
EXT FUNC 
FORM PAP 
REAL VAR 
LABEL 
I NT 4 VAR 
LABEL 
EXT FUNC 


0000 ERRORS 


ORBINAL PAG' 
ipf; POOR 


C >- 


W 




